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ABSTRACT 



This report discusses the theory, design, implementation and testing of a 
personal computer-based Multi-Frequency Modulation (MFM) packet commu- 
nications system. Transmitter,/ receiver programs provide software drivers for 
DA and A, D boards and perform symbol encoding, modulating, demodulating 
and decoding. The design and construction of a polarity coincidence correlator for 
receiver packet synchronization is presented. Experimental results show that the 
implemented MFM communication system conforms to theoretical analysis with 
acceptable bit error. Results also show that MFM can be uniquely adapted to a 
specific channel. 



THESIS DISCLAIMER 



The reader is cautioned that computer programs developed in this research 
may not have been exercised for all cases of interest. While every effort has been 
made, within the time available, to ensure that the programs are free of compu- 
tational and logic errors, they cannot be considered validated. Any application 
of these programs without additional verification is at the risk of the user. 
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I. INTRODUCTION 



Reliability and flexibility are two fundamental advantages of digital commu- 
nications. They are achieved through the speed and processing power of digital 
integrated circuits. Multi-Frequency Modulation (MFM) uniquely harnesses this 
power found in the modern personal computer to encode, modulate, demodulate 
and decode digitally formatted information, for sample rates currently up to 
100,000 samples per second. With computer-to-computer communication links 
implemented in a variety of mediums, such as wire, optical fiber, radio frequency, 
or acoustic, MFM readily adapts to a given medium and can emulate most ex- 
isting signal modulation formats. 

The focus of this thesis is the implementation of a MFM packet communi- 
cations system using industry standard personal computers (PC). Because of the 
limited processing speed of the PC, modulation and demodulation of the packet 
is performed in longer than real time. During the course of this project, trans- 
mitter and receiver software were developed to process data for establishing a link 
between PCs using D A and A D data acquisition boards. A correlator was de- 
signed and constructed to synchronize the receiver to the transmitted packets us- 
ing a synchronization baud at the beginning of the packets. 

Analysis of the system phase response led to signal modifications which re- 
sulted in improved signal quality. Signal-to-noise performance was evaluated on 
the system for various packet constructions. The scope of the report encompasses 
the theory of MFM. a description of its specific development and implementa- 
tion, and a summary of the evaluation of the system's performance. 
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II. THEORY OF MULTI-FREQUENCY MODULATION 



A. PACKET CONSTRUCTION 

The MFM signal set consists of "packets" of multiple tones which are ampli- 
tude and/or phase modulated. These tones are present simultaneously during a 
subinterval of the packet known as a baud. Packets can be located arbitrarily in 
the frequency spectrum and time as seen in Figure 1. 

The follow definitions are used in MFM [Ref. 1: pp. 5-6]: 



T 

AT 

k x 

L 

At 

/,= 1/Af 
A/= 1 jAT 
K 



Packet length in seconds 

Baud length in seconds 

Baud length in number of samples 

Symbol set. f, k is the phase of the k ,h tone in the l Th baud 

Number of baud per packet 

Time between samples in seconds 

Sampling or clock frequency for D/A and A/D conversion in Hz 
Frequency spacing between tones 
Number of MFM tones 



To ensure the packet can be uniquely represented by the k x samples, the 
Sampling Theorem states that f x , the sampling frequency, must be greater then 
twice the highest frequency in the signal set [Ref. 2: pp. 46-56]. Since f x = k x Af 
and Af is fixed by AT, the highest tone is kj 2 -1. Therefore, the signal set can 
consist of an arbitrary selection of tones from dc to fj 2 — Af. 

The analog representation for the MFM signal during the I th baud is 



k x !2 

Xj(t) = S )' i A lk cos{f.nkAft + f lk ), {l — \)AT < t < IAT. (1) 

k=0 



The first baud begins at t = 0 and the last baud, L, ends at T = LAT . Sampling 
x,{t) at t = nAt, where n is discrete time, and substituting A/=1/AT and 
At = AT/k x , produces a sampled output sequence of 
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Figure 1. MFM Signal Packet (after Ref. 1: p. 3.) 



k x l 2 

•*/(«) = Yj A ik cos ( + faX 0<,n^k x , (2) 

*=0 \ '* > 

a digital signal sequence of length k x samples. The k x point Discrete Fourier 
Transform (DFT) of is 



k x l2 

X ,(*') = £ Y k x A, t {e l **6{k' -k) + e- Jilt 6(k' - (k, - *))}, 0 <; *' < -1 .(3) 

*=o z 

From (3), it can be seen that the upper half of X \{k') is the complex conjugate 
image of the lower half, which is consistent with the symmetry property of a real 
sequence [Ref. 2: p. 402]. 
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B. SIGNAL GENERATION AND DEMODULATION 



In the past, generation of was attempted by constructing multiple phase 
lock loops, which were harmonically related. In December 1985, LCDR Deborah 
E. DeFrank, then at NPS, designed a coherent multifrequency synthesizer as a 
First step in producing [Ref. 3]. This method proved to be difficult to im- 
plement due to the phase jitter in the phase lock loops and changes required in 
hardware in order to change the harmonic frequencies used in the tone set. 

However, using the properties of DFT's, a baud of x,(n) can be generated by 
the host transmit computer by loading the First half of a complex valued array 
(0 to kJ2 —1) with the magnitude and phase of the tones to be included. To en- 
sure x,{n ) is real, the upper half, ( k x /2 ) to {k x — 1), is loaded with the complex 
conjugate of the values in the first half of the array at the image harmonics. An 
example of the DFT of a real x(n) having a period of k x = 16 and 3 tones is shown 
in Figure 2. The Inverse Discrete Fourier Transform (IDFT) generates the real 
discrete sequence, x,(n), which is clocked out of the computer thru a D, A con- 
verter, at f x samples per second. A signal packet is generated by L repetitions of 
the above process. 

Demodulation of MFM is the inverse of the signal generation process. The 
analog signal, x(r), is sampled at f x samples per second and converted to digital 
format with an A/D converter. The sampled values are loaded into the real 
components of a k x point complex array, with the imaginary component set to 
zero. A DFT is computed of the array to obtain the complex frequency repre- 
sentation of the sampled input. Since the upper half of the DFT is redundant 
information, only the lower half is retained for further processing. 

C. PROPERTIES OF MFM 

Several important properties of MFM signals are presented in this section. 



(I ) Orthogonality of signal during a baud. In continuous time 
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Figure 2. Data Structure (after Ref. 1: p. 9.) 
and in discrete time 



n= k x -1 

Yj X k( n M 0 = 

n = 0 



\m4k, 

0 



k = i 
k^i 



( 5 ) 



The derivation of the orthogonality of harmonically related signals is well known 
[Ref. 4: pp. 152-154]. Two advantages of an orthogonal signal set are: l)the 
noise, assuming additive white Gaussian noise(AWGN), affects each transmitted 
tone independently, thus simplifying statistical computations, and 2)no tone in a 
baud interferes with any of the other tones. This is shown by representing each 
tone in (1), neglecting phase and magnitude, as 



j*(f) = cos(2nkAft) rect(//A7 ) 



where 



rect(//A7) = 



{: 



Q^t^AT 

otherwise 



(6) 



V) 
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The Fourier transform of s k (t) is 

S k (J) = ATsmcff- kAJ)AT. (8) 

The spectrum of three adjacent tones, plotted in Figure 3, shows the peak of each 
tone coinciding with zero crossings of the spectrum of all other, thus minimizing 
their mutual interference. 

(2) Autocorrelation function (acf). Assuming x(n) is periodic, 
the circular acf of (2) is 



k x - 1 

r x(p)= ^ x{n)x{n@p), 0 < p < k x —l (9) 

n=0 

where © is a left circular shift. For a white bandlimited sequence with an even 
number of harmonics, the acf is given by 

■> sininKplkA 

r>) = (l/2).Vos(2, W si njig ( 1( » 

where k 0 = (k x + k 2 )l 2, the midband harmonic, and K = k { — k 2 + 1, the number 
of tones in the baud. Figure 4 shows the acf of a bandpass sequence with 
k x = 256. k x = 68 and k 2 = 83(A r = 16, k 0 = 75.5). These signal parameters con- 
form to those of the synchronization baud, which will be discussed later. Note the 
peak of the acf occurs at /?(0), this feature is fundamental to synchronizing the 
receiver to the incoming packet. 

(3) Matched filter. It is known that a matched Filter maximizes 
the signal-to-noise ratio for additive white noise [Ref. 4: pp. 88-89]. The DFT 
of ,x(«) at the frequency k is identical to the output of a filter matched to .x*(/r). 
This is due to the orthogonality property of the MFM signal and the linearity of 
the DFT [Ref. 1: pp. 17-21]. 

D. MODULATION 

Modulation is the process of encoding the source information onto a bandpass signal 
with a carrier frequency/. This bandpass signal is called the modulated signal 5(r), 
and the baseband source signal is called the modulating signal m(t) [Ref. 5: p. 204]. 
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CM 




The MFM signal can uniquely accomplish modulation in a number of ways. 
The signal ^/(«), as defined in (2), can be modulated in amplitude, frequency, and 
phase, by translating the message into changes in A, h ,k and <f), k respectively. Note, 
any combination of the modulation types is also possible, such as changing am- 
plitude and phase to produce quadrature amplitude modulation (QAM). The 
signal, 5(r), is called a bandpass signal. However, MFM can be bandpass or 
baseband and through multiplication with a carrier frequency .x,(rc) can be trans- 
lated to any frequency band desired. The signal sets considered in this thesis are 
bandpass and modulated using quadrature phase shift keying (QPSK) and dif- 
ferential quadrature phase shift keying (DQPSK). 
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Figure 4. ACF of a white bandpass sequence (after Ref. 1: p. 16.) 



( 1 ) QPSK. Conventional QPSK converts a digital input into 
four modulation voltage levels(symbols) to determine the phase of the transmitter 
output. To minimize the probability of symbol error the phases are spaced at 
multiples of rc/2 . A plot of the complex envelope of one tone of x,(«) is shown in 
Figure 5. The angle </> can take on the values of ± n/A and ± 3tt/4. 

Encoding MFM with QPSK is accomplished by loading the 
complex frequency domain array with the appropriate phase information. For 
example, if the digital input is '0110../. The first symbol, '01', would be loaded 
into the frequency bin, k, as Re[A f (/c)] = -AjJl and Im[A'(&)] = AjJl, where 
A k is the amplitude of tone k. Likewise, the second symbol TO' would be loaded 
into the second bin, k+ 1, as Re[A(/c + 1 )] = — A kArX lJl and 
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Imaginary 

(quadrature) 




Real 

(in phase) 



Figure 5. QPSK signal constellation. 

Im[A'(/c + 1)] = A k ^ljl. This would continue until all the tones in the baud were 
filled or there were no more symbols. As mentioned earlier, the discrete time do- 
main signal is produced by taking the IDFT of the complex frequency array. 
Each successive baud is encoded similarly. 

Decoding QPSK directly into bits is accomplished easily as fol- 
lows. Assuming a coherent receiver, decoding requires evaluating the polarity of 
the real and imaginary components of each frequency bin. Notice in Figure 5 
that the symbol mapping uses Gray encoding. This reduces the probability of 
bit error because errors caused by AWGN are likely to occur when the adjacent 
symbol is selected for the transmitted symbol; thus, the symbol error will contain 
only one bit error. Gray encoding also allows decoding straight into bits, with the 
right bit determined by the polarity of the real component and the left bit by the 
polarity of the imaginary component. The digital signal is obtained by succes- 
sively decoding each tone of each baud. 

(2) DQPSK. QPSK requires strict phase coherent regeneration 
of the sampling frequency to eliminate phase ambiguity. This results in a complex 
synchronization design dr distribution of the sampling clock frequency to each 
receiver. DQPSK provides a practical solution to the phase uncertainty problem 
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Previous Symbol Previous Symbol 






but at a cost of approximately 2.3 dB more in SNR in order to obtain the same 
probability of bit error as coherent QPSK [Ref. 6: p. 442]. Encoding MFM with 
DQPSK is similar to QPSK; however, DQPSK translates the original symbol set 
into a second "differential" symbol set, which is then encoded as QPSK. Trans- 
lated symbols are generated based on the input symbol and the previous trans- 
lated symbol. Figure 6 shows this translation. Notice that, regardless of the 
previous symbol, an input of '00' generates a new symbol in the same quadrant 
as the previous symbol. 'An input of '01' rotates the new symbol + n/2 radians 
from the previous, 'IF rotates n radians, and TO' rotates — n/2. In the receiver, 
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decoding is performed by determining the phase difference between successive 
pairs of tones, using complex arithmetic. The phase difference, A <£, is found from 

gf&t = 1)* _ g/A- L e (/A -M-i) 5 (11) 

and in the absence of noise will be 0, ± nj 2, or n radians. To realign the complex 
signal to the original constellations, is rotated by + rc/4 radians and the in- 
phase and quadrature components are decoded as with QPSK. 
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III. SYSTEM DEVELOPMENT 



Multi-Frequency Modulation joins the basic methods of digital communi- 
cations theory and signal processing. The theory behind MFM is not new [Ref. 
7], however, the ability to implement it inexpensively is linked to the recent de- 
velopment of inexpensive DSP chips for personal computers. The advent of 
packet switching in data communications also makes MFM a preferred choice for 
the MODEM because of its "packet" like format. 

A. SYSTEM DESCRIPTION 

An MFM communications system is shown in Figure 7. Information, de- 
noted by m{n), is modulated into a frequency band that will propagate over the 
available channel. The receiver converts the noise-corrupted signal, r(r), to an 
estimate of the source information, m(n). Theoretically, m(n) is any signal that 
can be represented in a digital format. To be practically useful, the conversion 
process must be real time. However, due to signal conversion speeds available for 
this research project, the input information had to be restricted to stored data 
Files. The ideal transmitter output, x(t), is described for the l' h baud by (1). 



n(t) 




Figure 7. MFM Communication System. 



B. BLOCK DESCRIPTION 
1. Transmitter 

The transmitter is subdivided into functional blocks, shown in Figure 8. 
As mentioned above, m(n) has been digitized and stored in a file. Therefore, the 
input is a string of binary digits. Processing is performed on a baud-by-baud 
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Figure 8. Transmitter Functional Block Diagram. 



basis, until the end of the data file, or maximum packet length has been obtained, 
whichever comes first. 

The encoder converts input symbols into complex values stored in the 
frequency domain array. The value of the symbols depends on the type of mod- 
ulation. For example, QPSK symbols are two bits long and are encoded as pre- 
viously discussed. The discrete signal, produced by computing the IDFT of the 
complex-valued frequency domain array, is loaded into the packet storage area. 
The controller determines the parameters of the packet based on the modulation 
type, baud, and message size; it then sequences the input data through the 
transmitter one baud at a time. Once the message has been processed, the entire 
stored digital signal, x(«), is transferred out at the selected rate, f x samples per 
second, through a D/'A converter. Depending on the channel's frequency re- 
sponse, filtering of the output may be desired to remove the high frequency am- 
plitude discontinuities introduced by the D/A converter. 

2. Receiver 

The receiver, shown in Figure 9, demodulates the MFM signal by re- 
versing the transmitter process. Data acquisition is the process that samples the 
analog signal at f x samples per second and converts it to digital format with an 
A/D converter. Though not shown in Figure 9, filtering of the input is recom- 
mended to bandlimit input noise and to reject out-of-band interference. The 
converted data sequence is stored until all data is acquired. Again, this is due to 
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time limitations in the signal processing algorithms available in the system used 
in this research task. The stored real values are accessed one baud at a time to 
perform a k x point DFT. The first half of the resulting complex values are de- 
coded to obtain the amplitude and phase modulation information. As in the 
transmitter, the controller sequences the data through the system. 




Figure 9. Receiver Functional Block Diagram. 

Notice the delay from input to output is the time it takes to acquire the 
entire packet and process the first baud. In contrast, real time processing must 
complete the k x point DFT in the same or less time than it takes to fill up one 
baud length buffer. Processing can alternate between buffers, and data will flow 
continuously through the system at the sampling rate. The storage requirements 
and processing time using f x = 61440 Hz, L =30, and k x =1024 samples are given 
in Table 1 for a system with a real time and non-real time DFT capability and 
for the non-real time system used in this research. The AT implementation 
processing time highlights the need for a high speed vector processor to perform 
the DFT computations. 

3. Synchronization 

In digital communications, various degrees of synchronization are re- 
quired. Typical degrees are: 

• Carrier synchronization, 

• Bit or symbol synchronization, 

• Word or frame synchronization. 
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Table 1. STORAGE REQUIREMENTS AND PROCESSING TIME 





Real time 


Non real time 


AT implementa- 
tion 


Storage required 
(# of samples) 


2048 


30720 


30720 


Time required for 
processing 1024 
samples (msec) 


< 16.7 


> 16.7 


10 x 10 3 



These have slightly different meanings and forms depending on the system. For 
the MFM system, carrier and frame synchronization may be required. 

Carrier synchronization is required if the modulation scheme is coherent, 
such as QPSK. This means generating an f x at the receiver in frequency and phase 
coherence with the transmitted f x . Though the packet does not contain the f x 
harmonic, every tone is harmonically related to it. The packet can be constructed 
to have a pilot tone separated from the modulation tones from which f x can be 
derived. For example, the modulation tones could range from k, to k M , giving a 
band of j + 1 consecutive tones. The pilot tone, k p , would then be placed several 
tones away to ensure minimum band interference from the message during its 
extraction. 

As mentioned previously, MFM, as implemented, is in the form of a 
transmission packet. To acquire the packet the receiver must know when to start 
sampling. This is accomplished by frame synchronization. Typically, unique 
words are inserted to mark the start of each frame. In MFM, the unique word is 
called a "synchronization baud", and it is added at the beginning of each packet. 
This baud is generated similarly to other bauds, except the tones and phases are 
predetermined. Acquisition of the received signal starts after successful detection 
of the synchronization baud. [Ref. 5: pp. 51 1-512, 8: pp. 293-295] 
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IV. SYSTEM IMPLEMENTATION 



In the previous chapters the tools and framework for developing MFM 
packet signals were established. One actual realization of this form of communi- 
cations will now be presented. Though the IBM Personal Computer (PC) was 
used, MFM can be implemented in a variety of ways. The PC, however, is ide- 
ally suited for MFM implementation for the following reasons: 

• Increased I/O channel maximum throughput rate of approximately 
100 KHz using direct memory access (DMA). 

• Digital signal processing allows encoding, decoding, modulation, demodu- 
lation, and channel equalization. 

• Signal processing algorithms are available in high level languages. 

• Packet construction is easily modified to conform to various channels. 

• External hardware is easily interfaced. 

• Cost is low. 

Successful data acquisition and decoding of (1) was the First and most im- 
portant goal of this thesis. Other goals affecting software and hardware develop- 
ment were the following: 

• Maximize packet size. 

• Develop synchronization circuitry. 

• Transmit, receive ASCII Files. 

• Develop software for testability and flexibility. 

• Develop software for statistical testing. 

A. SIGNAL PARAMETERS 

Software and hardware implementations can be easily modified to receive any 
packet construction. However, all software conforms to the signal parameters in 
Appendix A, established between NPS and NOSC for an acoustic application. 
This is a bandpass signal in the band from 16-20 KHz using a clock rate of 
61440 Hz. All baud lengths are powers of two, allowing utilization of 
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hardware software Fast Fourier Transform (FFT) algorithms to speed frequency 
transformation processing. 

Various baud sizes are available to give greater flexibility in adapting the 
packet to test specific channel parameters. In DQPSK, where the information is 
represented by the phase difference between adjacent tones, channel phase dis- 
tortion affects the shorter baud more, due to their larger A/. The maximum 
tone-to-tone phase error introduced by a linear phase channel is shown in 
Table 2. 



Table 2. PHASE ERROR FROM A LINEAR CHANNEL 





256 


512 


1024 


2048 


4096 


A/ (Hz) 


240 


120 


60 


30 


15 


Phase error 
(degrees) 


25.6 


12.8 


6.4 


3.2 


1.6 



Obviously the longer baud would be preferred when differential coding is between 
adjacent tones. However, if the channel introduces time-related distortion, like 
noise bursts, electrical glitches, or propagation fluctuations, a short duration baud 
is desired, and differential coding should be between the same tones on adjacent 
bauds. 

B. SOFTWARE 
1. Transmitter 

Generation of the MFM signal has been accomplished in previous re- 
search at NPS. Utilizing previous hardware and core software, this project ex- 
panded the transmitter software to provide encoding of 16-QAM and QPSK, 
encoded data files using DQPSK, constructed maximum size packets, generated 
a synchronization baud, and provided greater flexibility in initializing the DMA. 
Actual transmitter hardware will not be discussed; it has been the subject of 
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previous research [Ref. 9]. Two new transmitter programs, TRANSMIT and 
XMITMES, will be discussed. 
a. TRANSMIT 

The program TRANSMIT provides maximum flexibility in con- 
structing and encoding signal packets. It is an excellent tool for basic system 
testing and can be used as a training aid to demonstrate various digital modu- 
lation schemes. A procedural flow diagram of TRANSMIT is shown in 
Figure 10. SelectBaud is a user interactive procedure for establishing the pa- 
rameters of a 16-20 KHz packet and for choosing the encoding scheme to be 
utilized. Parameters k x and L are selected. From k x , the lower and upper 
bandlimits, k x and k 2 , are set. This generates a bandlimited packet with 
k 2 — k x + 1 tones. The data storage requirement is k x L, representing the total 
number of samples in the packet which are clocked out. The encoding scheme 
determines the path for encoding the complex frequency domain array. 




Figure 10. TRANSMIT algorithm. 



EncodeQPSK begins by displaying a four symbol QPSK constellation 
used as a reference in selecting the symbols over the band. The symbols for the 
K tones in the baud can be selected in a variety ways: 

• Symbols for all tones are randomly selected from a random generator. 

• The symbol for each tone is selected by the user. 

• Individual tones may be removed from the band. 
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This last feature allows construction of a baud with an arbitrary number of tones 
within the band determined by k l and k 2 . Selected symbols for the band are 
loaded into the complex frequency domain array with their complex conjugate 
image frequency. EncodeQAM is functionally the same as EncodeQPSK, except 
the symbols conform to a 16-QAM constellation [Ref. 1: p.27]. 

To obtain the real discrete time-domain sequence of the encoded 
baud, ComplexFFT computes the inverse FFT. ComplexFFT consumes the ma- 
jority of the processing time in the program due to the complex arithmetic oper- 
ations required, thus restricting the overall throughput of the system. For real 
time processing, the FFT algorithm must be accomplished by a hardware signal 
processor. 

Each value in the time domain sequence is represented as a real data 
type, occupying six bytes of memory. ScaleData converts these values down to a 
one byte format acceptable to the D/A converter and places them into a packet 
storage buffer. EncodeQPSK, ComplexFFT, and ScaleData are executed for 
each baud, until all L baud have been processed. To transmit the packet out of 
the computer, DMAINJT transfers samples at f x samples per second over the 
DMA channel to D A converter [Ref. 9,10]. 
b. XMITMES 

To demonstrate the suitability of MFM for transferring information 
from a source to a sink, the program XMITMES was written to transmit an 
ASCII file encoded using DQPSK. Affixed to the beginning of the packet is a 
synchronization baud. As shown in Figure 11, XMITMES has an even simpler 
structure then TRANSMIT, because all tones in the band are encoded using 
DQPSK. 

The synchronization baud is a predetermined sequence generated by 
SyncBaud. This baud is constructed as are all other baud, except k x is fixed at 
256, and tones 68 to 83 are encoded with the same random symbol pattern re- 
gardless of the packet construction or input message. This synchronization se- 
quence occupies the first 256 values in the packet buffer and therefore is the first 
to be clocked out of the computer. 
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Figure 11. XMITMES algorithm. 



As messages can be of various sizes, TailorPacket sets the maximum 
number of baud required for encoding. This is determined by dividing the num- 
ber of characters in the message by the number of characters that can be en- 
coded. 

DiffEncode encodes the message file into the complex frequency do- 
main array. It reads one character at a time; then breaks the eight bit character 
into four 2-bit symbols. The symbols are DQPSK encoded and stored in the 
frequency array. Once encoded, processing and signal output by ComplexFFT, 
ScaleData, and DMAINIT are the same as in TRANSMIT. 

2. Receiver 

Like the transmitter, the receiver requires hardware to interface the com- 
puter with the channel and processing software to demodulate the MFM signal 
packet. The channel interface is a high-speed data acquisition board, model 
DASH-16F, manufactured by MetraByte Corporation [Ref. 11]. Also at the re- 
ceiver, synchronization circuitry is required to detect the beginning of the packet. 
The synchronization circuitry is discussed in Section C. Implementation of the 
processing software assumes synchronization and a receiver f x the same as at the 
transmitter. Two receiver programs, RECEIVE and RECMES will be discussed. 
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a. RECEIVE 

The program RECEIVE, shown in Figure 12, processes and displays 
only the first baud of a packet. It is an excellent tool to experiment with the 
DASH-16F's software drivers and it provides a quick indication of the overall 
system performance. Tw r o procedures are available for data input, GetData and 
AcquireData. GetData reads in the time domain sequence file generated in the 
transmitter, eliminating all data acquisition hardware in the transmitter and re- 
ceiver; thus ensuring perfect synchronization and accurate evaluation of all en- 
coding, decoding, and processing software. The computational SNR of the system 
is determined under this condition, as the only noise present is from round off 
and truncation errors introduced in processing. 




Figure 12. RECEIVE algorithm. 

Analog data acquisition is performed by the procedure Acquiredata. 
It initializes and controls the DASH-16F using procedures written by Quinn- 
Curtis [Ref. 12]. AcquireData allocates memory to store the sampled values 
transferred from the board using the DMA controller. Direct Memory Access is 
the only data transfer mode capable of transferring data to memory at the re- 
quired f x , without disruption by other interrupt processes in the computer. Other 
important initialization parameters are triggering source and the number of sam- 
ples to be collected.The A/D may be triggered from two sources, a programmable 
interval timer or an external trigger source. The programmable interval timer 
divides either a 1 MHz or 10 MHz clock to derive the sampling rate of the trig- 
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ger. Since this method cannot produce an arbitrary f x , external triggering is used. 
After initialization, conversions take place on the positive transition of every 
trigger until the specified number of samples have been acquired and transferred 
to memory. Collected data is then converted into a format acceptable for further 
processing. 

DiffDecode determines the encoded symbols by differentially decod- 
ing the complex frequency array transformed from the sampled data by 
ComplexFFT. The output of the program is the frequency spectrum plot in Fig- 
ure 13. The lower plot is the spectral response of the kJ2 tones in the baud; the 
upper plot represents only tones from k x to k 2 . The color on the original display 
indicates the phase quadrant of a given tone. Graphing the data in this fashion 
provides quick qualitative analysis of the frequencies and their phase informa- 
tion. 




Figure 13. Baud magnitude spectrum. 
b. RECMES 

RECMES, shown in Figure 14, demodulates the ASCII encoded 
transmission produced by XMITMES. It differs from RECEIVE in that it can 
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Figure 14. RECMES algorithm. 

process a multiple baud packet. The user interactive procedure, PacketSetUp, 
tailors processing to the expected receive packet, using the inputs k x and L. 
AcquireData samples and stores a memory segment of data regardless of the 
packet size, Lk x . However, once stored, ConvertData, ComplexFFT, and 
DiffDecode process only L baud of the data. DiffDecode combines four 
differentially decoded symbols into one byte, representing the ordinal number of 
an ASCII character. To reconstruct the message, the characters are transferred 
to text file MESSAGE.DAT until processing is complete. For convenience, 
ShowMessage displays the recovered message. 

C. HARDWARE 

Synchronization of DQPSK MFM is obtained from a hardware correlator 
that is external to the host receiver computer. The 128 point correlator, illus- 
trated in Figure 15, provides the data acquisition board with sampling triggers 
synchronized with respect to time of arrival of the packet. Using only the polarity 
of the analog input, it functions as a matched filter to the last half of the 256 
point synchronization baud. This type of correlator is referred to as a polarity 
coincidence correlator(PCC). 

The hard limiter used to obtain the polarity information in the analog input 
is a fast, high precision, high gain, operational amplifier. During positive and 
negative portions of the input the output is +5 Vdc and 0 Vdc respectively. 
This unipolar signal is synchronized to the receiver's f x as it is clocked through a 
128 point serial shift register. To minimize bit instability, due to the shift register 
input being asynchronous to f x , the hard limiter slew rate should be as large as 
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Sampling 
trigger 
to receiver 



Figure 15. Polarity coincidence correlator. 



possible. For example, during testing, .15 bits out of 128 were unstable using an 
MCI 747 general purpose operational amplifier with a slew rate of 0.7 V/>sec. 
This was reduced to 2 bits using the faster LM318 with a slew rate of 
70 Vjfj. sec. 

The Is and 0's buses store the time-reversed polarity sequence of the syn- 
chronization baud. The voltage on the buses represent the correlation between the 
stored sequence and the sequence in the shift register. As the sequences come into 
alignment, the l's bus voltage increases while the 0's decreases. This inverse 
property is combined by a differential amplifier to give the total correlation. 
When both sequences match, the l's bus voltage is +5 Vdc and the 0's is 0 Vdc, 
giving a maximum differential voltage of 5 Vdc into the threshold detector. The 
threshold detector generates a synchronization trigger on detection of the corre- 
lation peak. A latch is set, enabling the D flip-flop to pass sampling triggers for 
data acquisition at the packet's beginning as required for demodulation. 

During design and testing the correlator progressed from a 16 to 128 point 
shift register. With each shift register output bus connected through a resistor 
voltage divider network, progressive testing was necessary to ensure power sup- 
ply current ratings were not exceeded. Maximal-length sequences were used as 
test correlation sequences due to their unique two-valued autocorrelation func- 
tions that are easily determined by finite-field arithmetic [Ref. 8: pp. 368-375]. 
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V. PERFORMANCE EVALUATION 



Initial testing of the MFM system was concerned with quickly determining 
the quality of the demodulated signals. Figure 13 provided an excellent repre- 
sentation of the demodulated and decoded signal and was used throughout sys- 
tem development. This plot uniquely displays the frequency domain magnitude 
and phase. Since QPSK symbols are at least 90 degrees apart each color (on the 
original display) represents a specific quadrant in the complex frequency plot. For 
example, a green line indicates a tone whose phase is in the second quadrant and 
whose magnitude is represented by the length of the line. Although each color 
represents a broad decision region in phase, it is still clear a phase shift of at least 
180 degrees occurs across the 4 KHz band. This shift is primarily due to the linear 
phase Butterworth filter used to simulate the system channel, and it contributed 
to the decision to develop a DQPSK encoded signal. 

A. SYSTEM FREQUENCY RESPONSE 

After development of all hardware and software, testing revealed significant 
phase fluctuations across the frequency band. The system phase response, plotted 
in Figure 16, was determined by subtracting the transmitted phase from the re- 
ceived phase for each tone. Ideally the response should be smooth. With 
DQPSK encoded from tone-to-tone, phase differences introduced by the system 
will result in reduced tolerance to additive noise. For example, notice the phase 
difference of 0.4 radians between tones 296 and 297. When noise is added, de- 
coding the symbol in error will occur at a lower noise level than it would if the 
phase difference was zero; thus a higher SNR is required for low error rate de- 
coding when the channel introduces tone-to-tone phase fluctuations. The tone- 
to-tone phase fluctuations were substantially reduced, as shown in Figure 17, by 
decreasing the magnitude of the encoded tones in the transmitter which moved 
the operating range of the D/A converter to a more linear region. Further anal- 
ysis of the effect of the D/A and A/D converters on phase shift is recommended 
in order to be able to specify them properly in future designs. 
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HRRMONIC NO. ( AF-60HZ ) 



Figure 16. Initial system response. 

A second important feature of the system phase response is its average slope 
across the band. The phase slope is affected by synchronization timing as shown 
in Figure 18. When the sequence through the correlator was delayed 2 bits, or 
2/61440 = 32.55 fisecs, relative to the analog input to the receiver, thus delaying 
data acquisition accordingly, the phase slope changes as predicted by the time 
delay Fourier transform theorem; 

x(t-T d ) X{f)e~ JajTd . (12) 

This shows the phase slope in radians for a signal delayed by T d seconds, is 2nT d 
radians/'Hz. For ’example, at /= Afk , , = (60)(300) = 1 8000 Hz and 
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HRRMONIC NO. ( AF-60HZ ) 



Figure 17. Improved system response. 

T d = 32.55 usees , the phase shift is 3.7 radians. This compares closely to the ac- 
tual phase shift of tone 300 in Figure 17 and Figure 18. 

In summary, system frequency response analysis identified phase shifts in- 
troduced by the system hardware and by the synchronization timing. The mag- 
nitude of the encoded tones and the synchronization sequence delay were selected 
experimentally in order to minimize tone-to-tone phase fluctuations and provide 
the comparatively flat response shown in Figure 18 for all subsequent testing. 
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HARMONIC NO. ( AF-60HZ ) 

Figure 18. System response with 2-bit synchronization delay. 

B. SNR PERFORMANCE 

The performance of the system when corrupted by additive noise is deter- 
mined by its output signal-to-noise ratio (SNR). For QPSK-encoded MFM the 
mean of the real and imaginary parts of the 2 K coefficients of the DFT represent 
the received signal amplitude, and their variance represents the noise power. The 
output SNR is defined as the ratio of the square of the mean to the variance of 
each for these 2 K coefficients [Ref. 1: p. 25]. 

Figure 19 shows the system SNR, which is the output SNR when there is no 
additive noise in the channel, versus the frequency spacing, A/, of the five baud 
types. As would be expected, DQPSK performs better with a smaller A/ because 
the phase difference between adjacent tones is smaller when A/ is smaller. 
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In Figure 20, the SNR out of the receiver is shown versus input SNR which 
was varied by adding noise to the analog signal. Theoretically, the output SNR 
should equal the narrowband input SNR in QPSK [Ref. 1: pp. 24-25]. The re- 
duction in output SNR at higher input SNRs is due to each baud approaching 
its maximum output SNR as set by system noise and shown in Figure 19. 
Table 3 lists the bit errors at various SNR input levels. There were approxi- 
mately 2500 bits transmitted for each entry. This clearly indicates the worsening 
effect of the channel on bauds with greater tone spacing A/ at high SNRs. 

A direct comparison of the MFM signal with commercial modems is difficult 
due to the wide variety of techniques that are used to achieve high-speed data 
transfer. However, with a bandwidth efficiency of 2 bits/s/Hz, a DQPSK MFM 
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Figure 20. SNR performance. 



signal can achieve a throughput of 6000 bits per second on a standard switched 
telephone line having a bandwidth of 3 KHz. Though not implemented in this 
project, several techniques could be utilized to increase the bandwidth efficiency, 
such as increasing the constellation size and/or using data compression algorithms 
as in commercial modems [Ref. 13, 14]. Using these techniques MFM could 
achieve a bandwidth efficiency of 4-8 Bits/'s/'Hz, which compares favorably to 
current high speed modems. 
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Table 3. BIT ERRORS IN 2500 BITS TRANSMITTED VS BAUD TYPE AND 



A/ 


k 


SNR (dB) 


K x 


5 


10 


15 


20 


240 


256 


165 


35 


18 


10 


120 


512 


154 


16 


3 


0 


60 


1024 


148 


17 


0 


0 


30 


204S 


140 


9 


0 


0 


15 


4096 


127 


10 


0 


0 
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VI. CONCLUSIONS AND FURTHER STUDY 



MFM is well suited to the signal processing environment of the personal 
computer. As implemented in this thesis, the packet signal provides a high speed, 
low bit error data transfer link between two industry standard computers. Syn- 
chronization hardware allows the link to use asynchronous data. 

The theory and properties of MFM have been discussed for a white bandpass 
packet signal. Packet construction and various modulation formats have been 
shown to be easily adapted to a given channel. The system was developed and 
implemented for testability and to be easily modified to accommodate a variety 
of applications. However, the data storage requirements and signal conversion 
speeds were unique to this project due to available hardware. As implemented, 
this system would require a multiple packet transmission to achieve practical op- 
eration. 

Analysis of the system's phase response led to adjustments in the encoded 
tone magnitudes and synchronization timing. This substantially improved the 
system's phase response. SNR test results indicate a smaller A/baud has superior 
performance in a linear phase channel. For Af less than 60 Hz, bit errors are 
acceptably low for input SNRs greater than approximately 15 dB. Exhaustive 
testing is required to compare bit error performance with conventional digital 
modulation schemes. 

Areas of further study should focus toward increasing the system bandwidth 
efficiency and implementation of hardware signal processing for real-time opera- 
tion. In fact, current research at NPS is directed at increasing the throughput 
speed of the system described. Improved bandwidth efficiency can be obtained 
through data compression algorithms and developing receiver software to. decode 
the 16-QAM encoded signal generated by TRANSMIT. Finally to provide a 
greater flexibility in adapting MFM to channels with significant uncompensated 
phase distortion, the option should be developed to encode DQPSK on a baud- 
to-baud basis. 
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APPENDIX A. DESIGN PARAMETERS 



Table 4. DESIGN PARAMETERS FOR A 1/15TH SECOND SIGNAL PACKET 
IN A 16-20KHZ BANDPASS CHANNEL. 



Baud length (sec) 


AT 


1/240 


1/120 


1/60 


1/30 


1/15 


No. of Bauds 


L 


16 


8 


4 


2 


l 


Tone spacing 


A/ 


240 


120 


60 


30 


15 


Lowest Harmonic 




68 


135 


269 


537 


1073 


Lowest tone freq 


fi 


16320 


16200 


16140 


16110 


16095 


Highest Harmonic 


h 


83 


166 


332 


664 


1328 


Highest tone freq 


h 


19920 


19920 


19920 


19920 


19920 


Samples per Baud 


k x 


256 


512 


1024 


2048 


4096 


Sampling freq 


fx 


61440 


61440 


61440 


61440 


61440 
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APPENDIX B. OPERATING INSTRUCTIONS 



A. PRELIMINARIES 
1. Hardware setup 
a. System setup 

Setup system as shown in Figure 21. All signal/trigger connections should 
be accompanied with a ground lead to minimize interference. 



r i bbon 
cab I e 




Figure 21. System interconnection diagram. 



b. D ASH- 16F switch setting: 



CHAN CNFG 8 

A/D - BIP 

DMA 1 
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GAIN (12.5v) dipswitch A, B--ON; C, D, USER-OFF 

Address (300h) dipswitch 9, 8-OFF; 7, 6, 5, 4-ON 



c. Trigger Requirements 

The D/A converter board uses negative logic as shown in Figure 22. The 
pulse width, PW, must be less then one nsec to ensure proper signal output. This is due 
to DMA request and acknowledge hardware. All hardware on the correlator and 
DASH-16F is positive edge-triggered. The trigger magnitude should be approximately 
5.3V peak, due to CMOS shift register chips in the correlator. 

transi i tter 



At= 1/6 1440 sec 
PU< Ipsec 



rece i ver 
PU — variable 
At= 1/6 1440 sec >1 

n 

Figure 22. Trigger specifications. 
d. Power Supplies 



+8Vdc supply HP6216A 30 v 500ma (or suitable substitute) 

-8Vdfc supply HP6216A 30v 500ma (or suitable substitute) 

+5Vdc supply HP6216A 30v 500ma (or suitable substitute) 



2. Software setup 
a. Transmitter 

1. Convert DMAIN1T.ASM and DMASTOP.ASM to BlNary files 
[Ref. 15: pp. 91-93]., 

2. Place files DMAINTT.BIN, DMASTOP.BIN, COMPFFT.INC, and 
FFT87B2.INC in the same directory as TRANSMIT. PAS AND XMITMES.PAS. 




+5.3U 

0 
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3. Compile to disk TRAXSMIT.PAS and XMITMES.PAS. 
b. Receiver 

1. Compile TP4D16.PAS to disk to create TP4D16.TPU. 

2. Place TP4D16.TPU in the same directory as RECEIVE.PAS and RECMES.PAS. 

3. Compile to disk RECEIVE.PAS and RECMES.PAS. 

B. MESSAGE TRANSMISSION PROCEDURE 

1. As necessary conduct preliminary equipment setup. 

2. Reset correlator (removes triggers to DASH-16F). 

3. Run XMITMES until "Ready to Transmit". 

4. Run RECMES (match baud size used in XMITMES). 

5. When receiver prompts "Ready to acquires", press enter key on transmitter. 

6. Demodulated message is displaved on receiver screen and stored in file 
MESSAGE.DAT. 

TRANSMIT and RECEIVE are similarly executed. However TRANSMIT would 
require manual loading of the synchronization baud in Table 5 by individual loading 
each tone with the appropriate symbol. 



Table 5. SYNCHRONIZATION BAUD SYMBOL SEQUENCE 



tone 


68 


69 


70 


71 


72 


73 


74 


75 


76 


77 


78 


79 


80 


81 


82 


83 


symbol 


3 


3 


2 


2 


4 


1 


2 


1 


2 


3 


4 


3 


1 


4 


3 


1 
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APPENDIX C. TESTING PROCEDURES 



A. RESPONSE TESTING PROCEDURE 

1. Conduct preliminary equipment setup. 

2. Run XMITMES and exit. 

3. Transfer file XMITDAT.DAT to the receiver computer. (Not required for subse- 
quent response testing if transmitted message and packet construction are un- 
changed.) 

4. Transmit message I AW message transmission procedure. (RECMES generates 
output file RECSTAT.DAT used by RESPONSE and STATS.) 

5. Run RESPONSE. (Ensure baud size is correctly edited into RESPONSE.) 

6. Graph output file RESP^an^.DAT. 

B. SNR TESTING PROCEDURE 

1. Conduct preliminary equipment setup. 

2. Add random noise source and true RMS voltmeter as shown in Figure 23. 



2 . 2KQ 




Figure 23. Test equipment interconnection. 



3. Set SNR level. 

4. Conduct steps 1-3 in the Response testing procedure. 

5. Run Statistics (results are in file STAT^au^.DAT). 
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APPENDIX D. CORRELATOR SCHEMATIC 
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Slates from each shift register 
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Latch and D Flip/Flop 



APPENDIX E. TRANSMIT 



program TRANSMIT; 

(-'Transmits variable length packets up to 61440. The length is 
determined by BCSTARRAY and can be take to a full segment of memory. 

No sync baud is used. All external files are assumed to be in the same 
directory*) 

type 

TNvector = array[ 0. . 4095] of real; 

TNvectorPtr = ATNvector; 

BCSTARRAY = array[ -28673. . 32767] of byte; (*stores output samples*) 



kx, (*baud size*) 

kl, (*lowest tone in band*) 

k2, (*highest tone in band*) 

NUMBAUDS, (-'number of bauds in packet*) 

BAUDCOUNT, (*current baud being processed*) 

BYTECOUNT, (-'number of bytes to be transferred by DMA*) 
I : integer; 

XREAL, (’’'complex freq/time domain arrays*) 

XIMAG : TNvectorPtr; 

QAM, (*encoding scheme*) 

INVERSE : boolean; (*direction of FFT*) 

ERROR : byte; (*status of ComplexFFT*) 

BCST : BCSTARRAY; (*packet storage buffer*) 

THEFILE : file of byte; (-'time sequence output file*) 

ANSWER :char; (*input variable*) 



(*$I FFT87B2. INC*) 
(*$I COMPFFT. INC*) 



(* *) 

procedure SelectBaud; 

(’•SelectBaud establishes QPSK or QAM encoding kx, kl, k2, and the number 
of baud in the packet*) 



var 

ANSWER : integer; 

MODANSWER : char; 

begin 

window(20, 10,80,20); 

(*Select encoding scheme*) 
repeat 
clrscr; 

writeln( 1 Select modulation desired 1 ); 
writeln; 

writelnC A --- QPSK’); 

writeln( ’ B 16 QAM’); 

read ln( MODANSWER); 
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until MODANSWER in f ' a' ' A' , ' b ' , ' B' 1 ; 
if MODANSWER in ['b\'B f j then 
QAM := true; 

( ,v Select buad size*) 
kx: =0; 
repeat 

clrscr; 

if kx < 0 then writeln('TRY AGAIN'); 

write ln( 'What is the length of the bauds (kx)?'); 

writeln( ' i. e. 256, 512, 1024, 2048, 4096'); 

readln( ANSWER); 

case ANSWER of 

256: kx: =256; 

512: kx: =512; 

1024: kx: =1024; 

2048: kx: =2048; 

4096: kx: =4096; 
end; (*case ANSWER*) 
if kx = 0 then kx := 
until kx > 0; 



-l; 



(*Set bandlimits*) 
case kx of 
256: begin 

kl: =68; k2: =83; 
end; 

512: begin 

kl: =135; k2: =166; 
end; 

1024: begin 

kl: =269; k2: =332; 
end; 

2048: begin 

kl: =537; k2: =664; 
end; 

4096: begin 

kl: =1073; k2: =1328; 
end; 

end; (*case kx*) 

(*Select number of baud*) 

write('How many bauds do your desire to transmit? 
readln( NUMBAUDS) ; 
window( 1,1,80,25); 
end; (*SelectBaud*) 

(■•- - — *) 

(* *) 

procedure DisplayQAM; 

(*DisplayQAM shows the 16-QAM constellation in the upper portion of the 
screen. A window is set at the bottom of the screen for further 
interaction with the user.*) 






begin 

clrscr; 
got oXY( 11,1); 
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write( chr( 179) ) ; 
write( chr( 179) ) ; 



write( chr( 218) ) ; 
for I:= 12 to 64 do 
write( chr( 196) ); 
write(chr( 191) ); 
for I:= 2 to 19 do 
begin 

GotoXY( 11,1); 

GotoXY( 65 , I) ; 
end; 

GotoXY( 11,20); 
write(chr( 192) ); 
for I := 12 to 64 do 
write(chr( 196)); 
write( chr( 217) ) ; 



window( 12,2,63, 19); ( 

write('This program encodes a 16 QAM multifrequency'); 



writeln( 

writeln( 



signal. ' ) ; 



writeln( ' 






' , chr( 179) , ' 






’) 


writeln( ' 


* 


* 


' ,chr( 179) , ' 


* 


* 


’) 

’) 


writeln( ' 


4 


3 


' , chr( 179) , ' 


2 


1 


writeln( ' 






' ,chr( 179) , ' 






') 

') 


writeln( ' 


* 


* 


' ,chr( 179) , ' 


* 


* 


writeln( ' 


8 


7 


' ,chr( 179) , ' 


6 


5 


’) 


writeln( ' 






' , chr( 179),' 






’ ) 


writeln( ' 
writeln( ' 


* 


* 


' , chr( 179) , ' 
' , chr( 179),' 


* 


* 


’) 

') 

') 


writeln( ' 


12 


11 


' , chr( 179) , ' 


10 


9 


writeln( ' 






’ ,chr( 179) , ' 






’) 


writeln( ' 


* 


* 


' , chr( 179) , ' 


* 


* 


') 

') 


writeln( ' 


16 


15 


' , chr( 179) , ' 


14 


13 


writeln( ' 






’ ,chr(179),' 






') 



(*active window*) 

window( 12,21,80,25); 
end; (*DisplayQAM*) 

(* *) 

(* *) 

procedure DisplayQPSK; 

(*DisplayQPSK shows the QPSK constellation in the upper portion of the 
screen. A window is set at the bottom of the screen for further 
interaction with the user.*) 



begin 

clrscr; 
gotoXY( 11,1); 
write( chr( 218) ) ; 
for I: = 12 to 64 do 
write( chr( 196) ) ; 
write(chr( 191) ); 
for I: = 2 to 19 do 
begin 

GotoXY( 11,1); write( chr( 179 )) ; 

GotoXY( 65 , I ); write( chr( 179) ); 
end; 
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GotoXYf 11,20); 
writef chrf 192) ); 
for I := 12 to 64 do 
writef chr( 196) ) ; 
writef chr( 217 ) ); 



windowf 12,2, 63,19); 
wr iteln( ' This program encodes 
writelnf 'The phase are shown 
writelnf ' 

write Inf 1 * 

writelnf ' 2 

writelnf ' 
writelnf ' 
writelnf ' 
writelnf ' 


a QPSK multifrequency signal. ' ); 
for one frequency. ' ); 

' , chrf 179) , ' 

' , chrf 179) , ' * 

' jchrf 179) , ' 1 

' ,chrf 179), ' 

' , chrf 179) , ' 

' , chrf 179) , ' 

' , chrf 179 ), ' 


writelnf ' 


' , chrf 179) , ' 


writelnf ' 


' ,chrf 179), ' 


writelnf ' 


' , chrf 179 ), ' 


writelnf 1 


' ,chrf 179) , ' 


writelnf' * 


', chrf 179),' 


writelnf' 3 


' , chrf 179) , ' 4 


writelnf ' 


' , chrf 179), ' 



') 

') 

') 

') 

') 

') 

') 

') 

') 

') 

') 

') 

') 

') 



(*active window*) 

windowf 12,21,80,25); 
end; (*DisplayQPSK*) 

C* *) 

(* *) 

procedure SelectQAM; 

(*SelectQAM selects and encodes the symbols for each tone in the band set 
by kl and k2. All symbols can be randomly or manually selected.*) 



RESPONSE : char; 

ANSWER, I : integer; 

VECTORARRAY : array[ 0. . 4095] of integer; 



C* *) 

procedure EncodeData; 

(*EncodeData loads the frequency domain arrays XREAL and XIMAG with the 
QAM symbols in VECTORARRAY and their complex conjugate image.*) 

var 

J : integer; 

TEMPR, TEMPI : real; 

begin 

f il lchar(XREALA , sizeof ( XREAL a ) , 0); 
f il lchar(XIMAG a , sizeof (XIMAG a ) , 0); 
for J: = kl to k2 do 
begin 

case VECTORARRAY[ J] of 
0 : TEMPI: = 0.0; 

1. . 4 : TEMPI: = 90. 0; 

5.. 8 : TEMPI: = 30. 0; 
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9. . 12 : TEMPI: =-30. 0; 

13. . 16: TEMPI: =-90. 0; 
end; (’-case VECTORARRAY*) 
case VECTORARRAY[ J] of 
0 : TEMPR: = 0. 0; 

1,5,9,13 : TEMPR: = 90. 0 
2,6, 10, 14: TEMPR: = 30. 0 
3,7,11,15: TEMPR: =-30. 0 
4,8, 12, 16: TEMPR: =-90. 0 
end; (*case VECTORARRAY*) 



(*Load complex conjugate image*) 
XREALa[ J] : = TEMPR; 
XREALa[ kx- J] : = TEMPR; 
XIMAGa[ J] : = TEMPI; 



XIMAGa[ kx-J] : = -TEMPI; 
end; (*J:=kl*) 
end; (*EncodeData*) 

(* 



begin 

repeat 

clrscr; 

writelnC' Select one of the following for baud 1 ,BAUDCOUNT, ' . ' ); 
writeln; 

writeln( 'R Randomly select all ' ,k2 - kl + 1,' vectors'); 

writeln(' I Individually select vectors'); 

readlnC RESPONSE); 

until RESPONSE in [ ' r ' , ' R ' , ' i ' , ' I ' ] ; 
f i 1 1 ch a r C VE CTOR ARR AY , s iz eo f C VECTOR ARR AY ) , 0 ) ; 



case RESPONSE of 

' r' , ' R ' : for I: = kl to k2 do 

VECTOR ARRAYf I] : =randomC 16)+1; 

' i ' , ' I ' : begin 

I: = kl; 

while I <= k2 do 
begin 

repeat 
clrscr; 
writelnC ' 
writelnC '# 
writelnC ' 17 
writelnC ' 18 
write( ' 19 
read In C ANSWER); 
until ANSWER in [ 1. . 19] ; 
case ANSWER of 



Vector Selection Menu'); 
vector desired'); 
random vector'); 
no vector' ); 

no more tones Tone ' ); 



1. . 16: VECTORARRAYf I] : = ANSWER; 

17: VECTORARRAY[ I] : = randomC 16)+1; 
18: ; C*nop*) 

19: I: = k2; 

end; C*case ANSWER*) 

I: =1+1; 

end; C*while I*) 
end; (*'i' , ' I '*) 
end; C*case RESPONSE*) 
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EncodeData; 
end; (*SelectQAM*) 

(" *) 

(* *) 

procedure SelectQPSK; 

(*SelectQPSK selects and encodes the symbols for each tone in the band 
set by kl and k2. All symbols can be randomly or manually selected.*) 



RESPONSE : char; 

ANSWER, I : integer; 

PHASEARRAY : array[ 0. . 4095] of integer; 



(* *) 

procedure EncodeData; 

(*EncodeData loads the frequency domain arrays XREAL and XIMAG with the 
QPSK symbols in VECTORARRAY and their complex conjugate image.*) 

var 

J : integer; 

TEMPR, TEMPI : real; 

begin 

fillchar(XREALA,sizeof(XREALA) ,0); 
f illchar(XIMAGA , sizeof (XIMAG a ) , 0); 
for J: = kl to k2 do 
begin 

case PHASEARRAY] J] of 
0 : TEMPI: =0.0; 

1.2 : TEMPI: = 80.0; 

3.4 : TEMPI: =-80. 0; 

end; (*case PHASEARRAY*) 
case PHASEARRAY] J] of 

0 : TEMPR: = 0.0; 

1.4 : TEMPR: = 80. 0; 

2.3 : TEMPR: =-80. 0; 
end; (*case PHASEARRAY*) 

XREALa] J] : = TEMPR; 

XREAL/ [kx-J] := TEMPR; 

XIMAGa] J] := TEMPI; 

XIMAGa] kx-J] := -TEMPI; 

end; (*J:=kl*) 
end; (*EncodeData*) 

c* *) 

begin 

repeat 

clrscr; 

writeln( ' Select one of the following for baud ' ,BAUDCOUNT, ' . ' ); 
writeln; 

writeln('R Randomly select all ' ,k2 - kl + 1,' phases'); 

writeln( ' I Individually select phases'); 

readln( RESPONSE); 

until RESPONSE in [ ' r ' , ' R ’ , ' i ’ , ' I ’ ] ; 
f i 1 1 cha r ( PHASEARRAY , s iz eo f ( PHASEARRAY ) , 0 ) ; 
case RESPONSE of 

' r ' , ' R' : for I: = kl to k2 do 
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Phase Selection Menu'); 
phase desired' ); 
random phase'); 
no phase' ); 

no more tones Tone ',1, 



PHASEARRAY[ I] : =random(4)+l; 

' i ' , ' I ' : begin 

I: = kl; 

while I <= k2 do 
begin 

repeat 

clrscr; 
writelnC ' 
writelnC ' # 
writelnC ' 5 
writelnC ' 6 
writeC ' 7 
readln C ANSWER); 
until ANSWER in [ 1. . 7] ; 
case ANSWER of 

1. . 4: PHASEARRAY[ I] : = ANSWER; 

5: PHASEARRAY[ I] : = randomC4) + l; 
6: ; (*nop*) 

7: I: = k2; 

end;C*case ANSWER*) 

I: =1+1; 

end; (*while I*) 
end; (*' i' ,'l'*) 
end;C*case RESPONSE*) 

EncodeData; 
end; C*SelectQPSK*) 

(* 



*) 



procedure ScaleData; 

(*ScaleData converts each real value in array XREAL down to a byte and 
stores the byte in packet storage buffer BCST. INDEX sets the location 
in the buffer of each byte*) 



INDEX, J, TEMP -.integer; 
DATA : byte; 



126); 



28673; 



begin 

for J : = 0 to kx-1 do 
begin 

TEMP := roundCXREALA[ J] 
if TEMP < 0 then 
TEMP := 0; 

DATA := TEMP; 

INDEX : = J + (BAUDCOUNT - 1) * kx 
BCST[ INDEX] : = DATA; 
writeCTHEFILE ,DATA); 
end;C*for J*) 
end; C*ScaleData*) 

(* - - *) 

c* — - *) 

procedure Dmainit(var BCST : BCSTARRAY; 

BYTECOUNT : integer); 

C*Assembly language procedure used to initialize and unmask the DMA for 
data transfer. The source must be converted to a BIN file.*) 
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procedure Dmastop; (*Masks DMA, stopping data transfer*) 
external ' DMASTOP. BIN 1 ; 

(* *) 

procedure TINPUT; 

(*This procedure reads in binary file, H0UTFA1.DAT, into an array, 
performs an FFT using ComplexFFT and displays the results 
with GraphData*) 

var 

NUMPTS : integer; 

C* *) 

procedure GetData; 

(*GetData reads time domain sequence from file H0UTFA1.DAT into XREAL*) 
var 

THEFILE : file of byte; 

DATA : byte; 

I, 

TEMP : integer; 

rdata : real; 

begin 

ass ign( THEFILE, 'H0UTFA1.DAT' ); 

reset(THEFILE); 

new (XREAL); 

new(XIMAG); 

f ill chart XREAL a , sizeof (XREALa ) , 0); 
f illchar(XIMAGA ,sizeof (XIMAGa) , 0); 
for I : = 0 to kx-1 do 
begin 

read( THEFILE , DATA); (*read file one byte at a time*) 

TEMP := DATA; 

temp := temp and 255; (*reduces number of bits used*) 

(*to represent the input data*) 

XREALa[I] := (TEMP - 126); (*loads XREAL array*) 
end; (*for I*) 
close(THEFILE); 
end; (*GetData*) 



C* - *) 

(* *) 

procedure PutData; 

(*Writes frequency domain arrays XREAL and XIMAG to file TFFTOUT. OUT*) 
var 

I : integer; 



OUTFILE .-text; 
begin 

as sign( OUTFILE, 'TFFTOUT. DAT' ); 

rewrite (OUTFILE); 

for I : = 0 to kx-1 do 

writeln( OUTFILE, I,' ' , XREALa [ I] , ' ' ,XIMAGa[ I] ); 
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close(OUTFILE); 
end; (*PutData*) 



procedure GraphData; 

(*GraphData plots the kx/2 points of complex array XREAL and XIMAG. 
Symbols are decoded and represented by the color of each tones 
spectral line.*) 

var 

I,J, (*loop indices*) 

KXDIV2, 

COLOR, 

X, XZ, 

Y, YZ : integer; 



begin 

GraphColorMode; 
palette( 2) ; 
Graphbackground( 0 ) ; 
draw(30, 10,30, 170,1); 
draw(30, 170, 300, 170,1); 
draw( 30, 155, 300, 155,1); 
draw( 30,80,300,80,1); 



(*y-axis*) 

(*x-axis for full spectrum*) 
(*x-axis for zoom spectrum*) 



(*every other dot equals one bin at 
draw( 166 , 170 , 166 , 174 , 1); (*x-axis scale marks*) 



256*) 



(*x-axis scale marks*) 
(*y-axis scale marks*) 
(*y-axis scale marks*) 



draw( 196, 170, 196, 174,1); 
draw(30, 113,26, 113,1); 
draw( 30,23,26,23,1); 

I: = 23; 

(*Grouping of 4 horizontal lines to identify background color 
used for one of the symbols*) 
repeat 

dr aw( 3 1,1, 300, I ,0); 
draw(31, 1+1,300, 1+1,1); 
draw( 31 ,1+2,300,1+2,2) 
draw( 31,1+3,300,1+3,3) 

I: =1+12; 
until I > 75; 

KXDIV2 := round( kx / 2); 



(*Decode symbols and assign color*) 
for I := 0 TO KXDIV2 do 
begin 

if (XREALa[I] >= 0) and (XIMAGa[ I] > 0) then 



COLOR 



if (XREALa[I] 



COLOR 
if (XREAL 
COLOR 
if (XREAL 
COLOR 



= 0 



< 0) and (XIMAGa[I] > 0) then 

< 0) and (XIMAGa[I] <= 0) then 
>= 0) and (XIMAGa[I] <= 0) then 



(*zoom spectrum*) 

if (I >= kl) and ( I <= k2) then 
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begin 

XZ := round( ( ( I * 4096.0) / kx) - 1040); 

YZ := 80 - round(0.5*(sqrt(sqr(XREAMI]) + sqr(XIMAGA[ I] )) )); 
draw(XZ, 80, X?,YZ, COLOR); 
end; 

(*full spectrum*) 

X := round( I /KXDIV2 * 256 +30); 

Y :=170-round(20*ln(sqrt(sqr(XREALA[ I] +1)+ 
sqr(XIMAGA[ I] +1) ))/ln( 10)); 
draw(X, 170, X,Y, COLOR); 



end; 

end; (*GraphData*) 

(* *) 

begin 

INVERSE := false; (*sets forward FFT*) 

ERROR := 0; 

GetData; 

NUMPTS := kx; 

writelnC ' Performing FFT 1 ); 

ComplexFFT ( NUMPTS , INVERSE ,XREAL , XIMAG , ERROR) ; 
writeln( ' Error = ERROR, ' hit the enter key ' ); readln; 

PutData; 

GraphData; 
end; (*Tinput*) 

(* *) 

(* *) 



begin 

repeat 

clrscr; 

as s ign( THEFI LE , ' H0UTFA 1 . DAT ' ) ; 
rewr ite(THEFILE); 
new(XREAL); 
new( XIMAG); 

INVERSE: =true; 

QAM; -false; 

SelectBaud; 

if QAM then 
DisplayQAM 
else 

DisplayQPSK; 

(*Packet contruction loop*) 

for baudcount : = 1 to numbauds do 
begin 

if QAM then 
SelectQAM 
else 

SelectQPSK; 

writelnC ' Performing I FFT 1 ); 
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Comp 1 exFFT( kx , INVERSE , XREAL , X I MAG , ERROR ) ; 
ScaleData; 

end; (*for BAUDCOUNT*) 

dispose(XREAL); 

dispose(XIMAG); 

close(THEFILE); 

BYTECOUNT := NUMBAUDS*kx-l; 

writeln( ' Press return to transmit' ); readln; 

Dmainit( BCST , BYTECOUNT) ; 

repeat 

writeln( ' Transmit some more? ( ,v yes or no*) 
readlnC ANSWER); 

until ANSWER in [ ' n' , * N' , 1 y' , ' Y' ] ; 
dmastop; 

window( 1,1,80,25); 
until ANSWER in ['n','N']; 

Tinput; 

end. 
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APPENDIX F. XMITMES 



program XMITMES; 

(*Transmits a syncbaud and message from file 'MESSAGE.DAT'. The message 
is encoded using QPSK. 'MESSAGE.DAT' is a text file. It should already 
exist before using this program. Output is used to collect data for 
TESTING*) 

const 

FIRST_ELEMENT = -28929; 
type 

TNvector = array[ 0. . 4095] of real; 

TNvectorPtr = ATNvector; 

(*Sync + 61440*) 

BCSTARRAY = array[ FIRST_ELEMENT. . 32767] of byte; 



kx, 

kl ,k2 , 1 ,w, 

NUMBAUDS ,MAXNUMBAUDS , 
BAUDCOUNT , BYTECOUNT , 
SYMBOLCOUNT , MAXNUMCHAR , 



MESSAGESIZE 
MAGNITUDE, 
CHARACTERS_PER_BAUD 
XREALjXIMAG 
INVERSE 

TEMPBYTE, ERROR 
BCST 
BYTEFILE 
TESTFILE 
ANSWER , 

NEXTCHAR 



: integer; 

: real; 

: TNvectorPtr; 

: boolean; 

: byte; 

: BCSTARRAY; 

: file of byte; 
: text; 

: char; 



(*$I FFT87B2. INC*) 
(*$I COMPFFT. INC*) 



c* - *) 

procedure SyncBaud; 

(*Process the synchronization baud and stores the 256 point time domain 
sequence at the beginning of the packet storage area. *) 



TNvectorPtr; 
integer; 
byte; 
real; 

begin 

new( SYNCREAL); 
new( SYNCIMAG); 

SYNC MAG: = MAGNITUDE; 



SYNCREAL, SYNCIMAG 
J, TEMP 
SYNCDATA 
SYNCMAG 
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(*load synchronization symbols*) 
f illcharC SYNCREAL a , sizeof ( SYNCREAL a) , 0); 
f illchar( SYNCIMAG a , sizeof ( SYNCIMAG a) , 0); 



SYNCREALa [ 68] : = -SYNCMAG; 
SYNCREALa [ 69] : = -SYNCMAG; 
SYNCREALa [ 70] : = -SYNCMAG; 
SYNCREALa [ 71] : = -SYNCMAG; 
SYNCREALa [ 72] : = SYNCMAG; 
SYNCREALa [ 73] : = SYNCMAG; 
SYNCREALa [ 74] : = -SYNCMAG; 
SYNCREALa [ 75] : = SYNCMAG; 
SYNCREALa [ 76] : = -SYNCMAG; 
SYNCREALa [ 77] : = -SYNCMAG; 
SYNCREALa [ 78] : = SYNCMAG; 
SYNCREALa [ 79] : = -SYNCMAG; 
SYNCREALa [ 80] : = SYNCMAG; 
SYNCREALa [ 8 1] : = SYNCMAG; 
SYNCREALa [ 82] : = -SYNCMAG; 
SYNCREALa [ 83] : = SYNCMAG; 



SYNCIMAGa[ 68] 


= 


-SYNCMAG 


SYNCIMAGa[ 69] 


= 


-SYNCMAG 


SYNCIMAGa[ 70] 


= 


SYNCMAG 


SYNCIMAGa[ 71] 


= 


SYNCMAG 


SYNCIMAGa! 72] 


= 


-SYNCMAG 


SYNCIMAGa[ 73] 


= 


SYNCMAG 


SYNC IMAG a [ 74] 


= 


SYNCMAG 


SYNCIMAGa[ 75] 


= 


SYNCMAG 


SYNCIMAGa [ 76 j 


= 


SYNCMAG 


SYNCIMAGa[ 77] 


= 


-SYNCMAG 


SYNCIMAGa[ 78] 


= 


-SYNCMAG 


SYNCIMAGa[ 79] 


= 


-SYNCMAG 


SYNCIMAGa[ 80] 


= 


SYNCMAG 


SYNCIMAGa[ 81] 


= 


-SYNCMAG 


SYNCIMAGa[ 82] 


= 


-SYNCMAG 


SYNCIMAGa[ 83] 


= 


SYNCMAG 



(*complex conjugate image*) 
for J := 68 to 83 do 
begin 

SYNCREALa[ 256-J] : = SYNCREALa [ J] ; 
SYNCIMAGa[ 256-J] : =-SYNCIMAGa[ J] ; 
end; (*for J*) 



Comp lexFFT( 256 , INVERSE , SYNCREAL , SYNC IMAG , ERROR) ; 

(*scale data/load time sequence*) 
for J : = 0 to 255 do 
begin 

TEMP: =round( SYNCREAL a [J] + 126); 
if TEMP < 0 then 
TEMP: =0; 

SYNCDATA: =TEMP; 

BCST[ J+FIRST_ELEMENT] : =SYNCDATA; 
end;(*for J*) 
dispose( SYNCREAL); 
dispose( SYNCIMAG); 
end; (*SyncBaud*) 

(* - *) 

(* *) 

procedure SelectBaud; 

(*SelectBaud establishes kx 3 kl, and k2*) 
var 

ANSWER : integer; 



begin 

kx: =0; 

(*select baud size*) 
repeat 

if kx < 0 then writeln('TRY AGAIN'); 

writeln( ' What is the length of the bauds (kx)?'); 
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writeln( ' i. e. 256, 512, 1024, 2048, 4096'); 
read ln( ANSWER); 
case ANSWER of 
256: kx: =256; 

512: kx: =512; 

1024: kx: =1024; 

2048: kx: =2048; 

4096: kx: =4096; 
end; (*case ANSWER*) 
if kx = 0 then kx := -1; 
until kx > 0; 

(*set tone limits*) 
case kx of 
256: begin 

kl: =68; k2: =83; 
end; 

512: begin 

kl: =135; k2: =166; 
end; 

1024: begin 

kl: =269; k2: =332; 
end; 

2048: begin 

kl: =537; k2: =664; 
end; 

4096: begin 

kl: =1073; k2:=1328; 
end; 

end; (*case kx*) 
end; (*SelectBaud*) 

C* 

(* - 

procedure TailorPacket; 

(*TailorPacket sets the maximum number of baud required to encode the 
message*) 

begin 

MESSAGESIZE: = f ilesize(BYTEFILE); 

writeln( ' Message is ', MESSAGESIZE , ' bytes,'); 

(*kx/2 is the number of bit/baud for QPSK. kx/2-2 is the number 
for DQPSK. Each character is 8 bits*) 

CHARACTERS_PER_BAUD: =(kx/8 - 2)/8; 

(*61440/kx is maximum number of bauds possible*) 

MAXNUMCHAR: = trunc( 61440. 0/kx * CHARACTERS_PER_BAUD ) ; 
if MESSAGESIZE > MAXNUMCHAR then 
begin 

writeln( ' Message is to large. The last ', 

MESSAGESIZE - MAXNUMCHAR, 

' characters will not be transmitted. '); 

MESSAGESIZE: ^MAXNUMCHAR; 
end; 



MAXNUMBAUDS: =trunc( MESSAGESIZE / CHARACTERS_PER_BAUD) ; 



(*ensure last few characters are included*) 
if frac( MESSAGES I ZE / CHARACTERS_PER_BAUD) > 0. 0 then 
MAXNUMBAUDS: =MAXNUMBAUDS + 1; 

repeat 

writeln; 

writeln( 'Enter number of ' ,kx, ' bauds to process. ' 
MAXNUMBAUDS,' is the maximum.’); 
readln( NUMBAUDS ) ; 

until NUMBAUDS in [ 1. . MAXNUMBAUDS] ; 
end; (*TailorPacket*) 



procedure DiffEncode; 

(*DiffEncode differential encodes symbols on a tone-to-tone basis. 
BYTEFILE is read from one byte at a time. The byte is isolated into 
2-bit groups and stored in BITS. BITS is then used to DQPSK encode the 
frequency domain arrays XREAL and XIMAG. Bytes partially encoded are 
carried over into the next baud by global variable TEMPBYTE. *) 

var 

J : integer; 

BITS rbyte; 

begin 

fillchar(XREALA,sizeof(XREALA) ,0); 
fillchar(XIMAGA,sizeof(XIMAGA),0); 

(*first tone of every baud set to pi/2*) 

XREALa[ kl] : = MAGNITUDE; 

XIMAGa[ kl] : = MAGNITUDE; 

if SYMBOLCOUNT = 0 then 

read(BYTEFILE, TEMPBYTE); 

(*break apart character*) 
for J: = (kl + 1) to k2 do 
begin 

SYMBOLCOUNT: =SYMB0LC0UNT + 1; 
if frac( SYMBOLCOUNT / 4) = 0.25 then 
BITS: = (TEMPBYTE and $C0) shr 6; 
if frac( SYMBOLCOUNT / 4) = 0. 5 then 
BITS: = (TEMPBYTE and $30) shr 4; 
if frac( SYMBOLCOUNT / 4) = 0.75 then 
BITS: = (TEMPBYTE and $0C) shr 2; 
if frac( SYMBOLCOUNT / 4) = 0. 0 then 
begin 

BITS: = TEMPBYTE and $03; 
if not E0F( BYTEFILE) then 
read( BYTEFILE, TEMPBYTE) 
else 

TEMPBYTE: =$40; (*fill character*) 

end; 

if (BITS < 0) and (BITS > 3) then 

writeln( ' Bits not assigned properly'); 
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(*dif f erential encode*) 
case BITS of 



0: begin XREALa [ 
XIMAGa 
1: begin XREALa 
XIMAGa 
2: begin XREALa 
XIMAGa [ 
3: begin XREALa [ 
XIMAGa [ 

end; (*case BITS*) 
end; (*for J*) 



J] := XREALa [ J-l] ; 
J] := XIMAGa [J-l] ; 
J] :=-XIMAGa[ J-l] ; 
J] := XREALa [J-l] ; 
J] := XIMAGa [J-l] ; 
J] :=-XREALa[ J-l] ; 
J] :=-XREALa[ J-l] ; 
J] :=-XIMAGa[ J-l] ; 



end; 

end; 



end; 

end; 



(*complex conjugate image*) 
for J: = kl to k2 do 
begin 

XREALa [kx - J]:= XREALa[J]; 

XIMAGa [ kx - J] :=-XIMAGa[ J] ; 

writeln(TESTFILE j BAUDCOUNT: 4, J: 5 , trunc( XREAL a [ J] ): 6, 
trunc(XIMAGA[ J] ): 6); 



end; 

end; (*Diff Encode*) 



procedure ScaleData; 

(*ScaleData converts each real value in XREAL down to a byte and stores 
the byte in the packet storage buffer, BCST. INDEX establishes the 
location in the buffer of each byte in the packet.*) 



INDEX, J, TEMP 
DATA 



: integer; 
: byte; 



begin 

for J : = 0 to kx-1 do 
begin 

TEMP := round( XREALa [J] + 126); 
if TEMP < 0 then 
TEMP := 0; 

DATA := TEMP; 



(*256 is added to INDEX to start message bauds 
after the sync baud*) 

INDEX := J+( BAUDCOUNT- l)*kx+FIRST_ELEMENT+256; 

BCST[ INDEX] : = DATA; 
end; (*for J*) 
end; (*ScaleData*) 

(* *) 

(* *) 

procedure Dmainit(var BCST : BCSTARRAY; 

BYTECOUNT : integer); 

(*Assembly language procedure used to initialize and unmask the DMA for 
data transfer. The source code must be converted to a BIN file.*) 



external 'DMAINIT5.BIN 1 ; 

(* *) 
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procedure Dmastop; 

(*Masks DMA, stopping data transfer.*) 
external ' DMASTOP. BIN' ; 

(* *) 

(* *) 

begin 

clrscr; 

(^contains hex values to be encoded and transmitted*) 
ass ign(BYTEFILE, 'MESSAGE. DAT' ); 
reset(BYTEFILE); 

(*0utput file of encoded symbols. Used for system testing*) 
a s s i gn ( TE S TF I LE , ' XM I TD AT . D AT ' ) ; 
rewrite(TESTFILE); 

INVERSE: =true; 
repeat 

writelnC ' Enter magnitude of tones. (must be > 0.0) '); 
r ead ln( MAGNITUDE ) ; 
until MAGNITUDE > 0. 0; 

writeln( ' Loading sync baud.'); 

SyncBaud; 

SelectBaud; 

TailorPacket; 

SYMB0LC0UNT: =0; 

TEMPBYTE: =$00; 

writelnC ' Number of bauds is ' ,numbauds); 

new(XREAL); 

new(XIMAG); 

for baudcount : = 1 to numbauds do 
begin 

Dif fEncode; 

writelnC 'Performing IFFT ' , BAUDCOUNT, ' 

NUMBAUDS - BAUDCOUNT , ' left'); 

Comp 1 exFFT( kx , I NVERSE , XREAL , XIMAG , ERROR) ; 

ScaleData; 

end; (*for BAUDCOUNT* ) 

dispose(XREAL); 
dispose(XIMAG) ; 
close(BYTEFILE); 
close(TESTFILE); 

BYTECOUNT := 256 + NUMBAUDS*kx - 1; 
repeat 

writeln( ' Press return to transmit '); readln; 

Dmainit(BCST, BYTECOUNT); 
repeat 

writelnC ' Transmit some more? (*yes or no*) ' ); 
read In (ANSWER); 

until ANSWER in [ ' n' , ' N' , ' y ' , ' Y' ] ; 
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dmastop; 

until ANSWER in [ 'n' , ' N ' ] ; 

(* reset(TESTFILE); 

while not EOF(TESTFILE) do 
begin 

while not EOLN(TESTFILE) do 
begin 

read(TESTFILE.NEXTCHAR); 

write(NEXTCHAR); 

end; 

readln(TESTFILE); 

writeln; 

end; 

close(TESTFILE); *) 
end. 
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APPENDIX G. DMAINIT 



codeseg segment 



public dmainit 
assume cs: codeseg 



; procedure DMAINIT ( BLKADDRESS : XMITPOINTER; 

; BYTECOUNT : INTEGER); 

; this procedure initializes dma channel 3 and sets the 
; parameters to output the array best by passing the address 
; of the array start on the stack. BYTECOUNT is the number 
;of bytes to transfer and is pushed on the stack by the 
; calling program 

dma equ 0 

dmapage equ 80h 

dmainit proc near 



les di,dword ptr[ bp+6] ; move add of best into es: di 



push bp 
mov bp,sp 



;use bp to address stack 



mov al,5bh 

out dma+ll,al 

out dma+12,al 

mov ax,es 

mov cl ,4 

rol ax, cl 

push ax 

and al,0fh 

out dmapage+2,al 

pop ax 

and al,0f0h 

add ax,di 

out dma+6,al 

mov al,ah 

out dma+6,al 

mov ax , [ bp+4] 

out dma+7,al 

mov al,ah 

out dma+7,al 

mov a 1,3 

out dma+10,al 

pop bp 

ret 6 



; get page offset 

; output waveform buffer start addr 



; output dma byte count 



; pop 6 bytes off stack for addr of best 



; unmask ch 3 to start 



; store in ch 3 dma page reg 



;save ax for dma start addr 



; reset first/last ff 

; calc high order 4 bits of buffer area 



; dma chan 3 single mode, read, autoinit 



dmainit endp 
codeseg ends 
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APPENDIX H. DMASTOP 



codes eg segment 

public dmastop 
assume cs: codeseg 



procedure DMASTOP; 

this procedure stops dma channel 3 



dmastop 



dmastop 

codeseg 



dma equ 0 

dmapage equ 80h 

proc near 

push bp 

mov al,7 ; mask ch 3 

out dma+10,al 

pop bp 

ret 

endp 

ends 



to stop 
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APPENDIX I. RECEIVE 



program RECEIVE; 

(*This program acquires an analog signal, stores the raw data in a memory 
buffer, converts the data to TYPE real, performs an FFT using COMPFFT, 
differential decodes between adjacent tones and displays the result in 
the frequency domain with the phase quadrant represented by color. 

Data input can also be read in from the time domain sequence file 
H0UTFA1.DAT generated by program TRANSMIT*) 

uses Graph, Crt, tp4dl6; 

(*$I-*) 

(*$R-*) 

const 

kx = 256; (*set baudlength*) 

kxml = 255; 

type 

TNvector = array[ 0. . kxml] 

TNvectorPtr = ATNvector; 



var 

INVERSE : boolean; 

XREAL, XIMAG : TNvectorPtr; 

ERROR : byte; 

NUMPTS , 

kl,k2 : integer; 

ANSWER : char; 

(*$I FFT87B2. INC*) 

(*$I COMPFFT. INC*) 

(* *) 

(* *) 



procedure AcquireData; 

(*AcquireData intializes Metrobyte DASH -16F data acquisition board. 
Using TTOOLS procedures D16_init and D16_ainm. Data transfer is 
controlled by the DMA controller and initialized by D16_ainm, and 
disabled by D16_dma_int_disable. TTOOLS procedures are external 
procedures included by 'uses’ tp4dl6.*) 

const max_buffer = 1000; 

var i: integer; 

rate: real; 

cnt_num, mode, cycle, trigger, 
base_adr, err_code, int_ level, dma_ level, 
board_num, chanlo, 

op_type, status, next_cnt , err_code_s : integer; 



of real; (*TYPE for real and imaginary 
data for FFT routing*) 

{-'Pointer for FFT data array 
which allows dynamic allocation 
of memory*) 
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dmaPointer: pointer; 
datavector: Ainteger; 

ad_data, chan_data: array[ 0. . max_buffer) of integer; 

begin 

clrscr; 

(*al locates memory on the heap and returns a pointer to the start 
of the buffer*) 

GetDMABuf fer(max_buf fer , dmaPointer , err_code); 

(*This statement assigns a generic pointer to a variable of a 

specific pointer type, i. e. Ainteger, so that the pointer can be 
passed to the dl6_ainm routine.*) 
datavector := dmaPointer; 

board_num := 0; int_level := 7; dma_level := 1; 
base_adr := $300; 

(•'initialize the driver*) 

dl6_init(board_num,base_adr , int_level ,dma_level ,err_code); 
chanlo : = 0; 

cycle: =0; (*0 - one sweep of the DMA 1 - autoinitialize*) 

trigger: =0; (*0 - external 1 - internal*) 
cnt_num: =kx; (*# of samples*) 

rate := 10000.0; (*used for internal trigger*) 
mode := 2; (*DMA mode*) 

writeln( ' Ready to acquire'); 

(*colects kx analog values using DMA and stores in a buffer*) 
dl6_ainm( board_num , chanlo , mode , cycle , trigger , cnt_num , rate , 
datavectorA , err_code); 

(*status indicates the progress of acquisition. When all 
samples have been acquired status=0*) 
status := 11; 

(*wait until all data acquired*) 
repeat 

dl6_dma_int_st atus ( board_num , op_type , status ,next_cnt , err_code_s ) ; 
until status = 0; 
w T riteln( ' Data received'); 

if err_code <> 0 then 

dl6_print_error( err_code) 
else 
begin 

writeln( ' Processing data'); 

(•'converts left justified data, returns the true binary value 
of the sampled data*) 

dl6_convert_data(2047,cnt_num, datavector a , ad_data[ 0] , 
chan_data[0] ,0,err_code); 

new(XREAL); 

new(XIMAG); 

f i 1 lchar ( XREAL a , s izeo f ( XRE ALa ) , 0 ) ; 
fillchar(XIMAGA,sizeof (XIMAGa),0); 
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for I:= 0 to kxml do 
begin 

XREALa[I] := ad_data[ I] /5; 
end; 



(*stop DMA operation if in autoinitialization*) 
d I6_dma_int_dis able( board_num , err_code ) ; 



(*frees memory allocated with GetDMABuf fer*) 
FreeDMABuf f er(max_buf f er , dmaPointer , err_code) ; 



end; (*Acquire*) 

(* 

(* 

procedure GetData; 

(*GetData loads XREAL with the real time domain sequence generated by 
program TRANSMIT. Only one baud is loaded*) 



THEFILE 

DATA 

I, 

TEMP 



: file of byte; 
: byte; 



: integer; 



begin 

ass ign( THEFILE, 'H0UTFA1. DAT' ); 

reset(THEFILE); 

new (XREAL); 

new(XIMAG); 

f illchar(XREAL/\ ,sizeof (XREALa), 0); 
f illchar(XIMAGA , sizeof (XIMAGa), 0); 
for I : = 0 to kxml do 
begin 

read(THEFILE,DATA); (*read file one byte at a time*) 

TEMP : = DATA; (*puts byte into integer variable*) 

XREALa[I] := (TEMP - 126); (*loads XREAL array*) 

end; (*for I*) 
close(THEFILE); 
end; (*GetData*) 



C* 



*) 



*) 

procedure DiffDecode; 

(*DiffDecode differentially decodes complex frequency domain arrays XREAL 
and XIMAG. Four decoded symbols are recombined into a byte and 
transferred to file BYTESOUT.DAT.*) 



SYMBOLCOUNT : intege 

TEMPRE AL , TEMP I MAG : real; 

BITS ,TEMPBYTE : byte; 

TEMPO HAR : char; 

OUTFILE : TEXT; 



begin 
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assign(OUTFILE, 'BYTESOUT. DAT' ); 
rewrite(OUTFILE); 

SYMBOLCOUNT: =0; 

TEMP BYTE: =0; 

for I := kl to (k2 - 1) do 
begin 

(*Complex multiply two adjacent tones, I and the complex conjugate of 
1+1. This will give the phase difference between the two tones. 

The answer is in rectangular notation*) 

TEMPREAL: =XREALa[ I] * XREALa[ 1+1] + XIMAGa[ I] * XIMAGaJ 1+1] ; 
TEMPIMAG: =XREALa[ I] * XIMAGa[ 1+1] - XREALa[ 1 + 1] * XIMAGa[ I] ; 

(*Complex multiply (TEMPREAL + j TEMPIMAGE) and (1+j). This rotates 
the differential vector 4-5 degrees. XREAL and XIMAG are used to 
store the results. This eliminate the original data*) 

XREALa[ I] :=( TEMPREAL - TEMPIMAG)/ 80; 

XIMAGa[ I] : =(TEMPREAL + TEMPIMAG)/80; 



if (XREALa[I] 
BITS: =$00; 
if (XREALa[I] 
BITS: =$01; 
if (XREALa [ I] 
BITS: =$03; 
if (XREALa[I] 
BITS: =$02; 
SYMBOLCOUNT : 



>= 0) and (XIMAGa[I] > 0) then 

< 0) and (XIMAGa[I] > 0) then 

< 0) and (XIMAGa[I] <= 0) then 
>= 0) and (XIMAGa[I] <= 0) then 

= SYMBOLCOUNT + 1; 



(*fill TEMPBYTE with four symbols*) 
if frac( SYMBOLCOUNT / 4) = 0.25 then 
TEMPBYTE := (BITS shl 6); 
if frac( SYMBOLCOUNT / 4) = 0. 5 then 

TEMPBYTE := (BITS shl 4) or TEMPBYTE; 
if frac( SYMBOLCOUNT / 4) = 0.75 then 
TEMPBYTE := (BITS shl 2) or TEMPBYTE; 
if (frac( SYMBOLCOUNT / 4) = 0.0) then 
begin 

TEMPBYTE := BITS or TEMPBYTE; 

TEMPCHAR := chr( TEMPBYTE); 

write(0UTFILE,TEMPCHAR); 

TEMPBYTE: =0; 
end; (*if frac*) 
end; (*for I*) 



TEMPCHAR := chr( TEMPBYTE); 

wr ite(0UTF I LE, TEMPCHAR); (*puts two 0 at end of each baud*) 
close(OUTFILE); 



XREALa[ k2] : =0; 

XIMAGa[ k2] : =0; 
end; (*Dif fDecode*) 

(* - 

(* 

procedure GraphData(XR, XI: TNVectorPtr); 

(*GraphData graphs complex arrays XR and XI. Two scales are used: Full 
scale, 0 to kx/2 tones and Zoom scale, kl to k2*) 



-*) 

■*) 



63 



(*loop indices*) 



grDr iver , grMode , ErrCode , 

Bj, 

kl ,k2 , 

MXDIV2, 

COLOR, 

X, XZ, 

Y, YZ : integer; 



begin 

kl := round(kx * 67. 0 / 256.0); 
k2 := round(kx * 83.0 / 256.0); 
grDriver : = EGA; 
grMode := EGAHi; 

initgraph( grDriver , grMode, ' c: TP4 GRAPHICS’); 

ErrCode := GraphResult; 
if ErrCode <> grOK then 
begin 

writeln( ' Graphics error: ', GraphErrorMsg(ErrCode) ); 

writeln( ' Program aborted... '); 

Halt(l); 

end; (*if ErrCode*) 

SetBKCo lor (black); 

SetColor(white); 



(*axis contruction*) 

line( 60,300,600,300); 
line(68, 140,600, 140); 
line (33 2, 300, 3 32, 3 08); 
line( 392,300,392,308); 
line (60, 300, 60, 308); 
line(572,300,572,308); 
line(96, 140,332,300); 
line(576, 140,392,300); 



(*x-axis for full spectrum*) 
(*x-axis for zoom spectrum*) 
(*x-axis scale marks*) 
(*x-axis scale marks*) 



MXDIV2 := round(kx / 2); 

for I := 0 to MXDIV2 do 
begin 

if (XRa[I] >= 0) and (XIa[ I] > 0) then 
SetColor( lightmagenta) ; 

if (XRa[I] < 0) and (XIa[I] > 0) then 
SetColor( lightgreen); 

if (XRa[I] < 0) and (XIa[I] <= 0) then 
SetColor( lightred); 

if (XRa[I] >= 0) and (XIa[I] <= 0) then 
SetColor( yellow); 

(*zoom spectrum*) 

if (I > kl) and (I <= k2) then 
begin 

XZ := round( 2*( ( I * 4096.0 / kx) - 1040)); 

YZ := 140 - round( sqrt( sqr( XRa [ I] ) + sqr( XI a [I] ) ) ) ; 

line( XZ , 140 ,XZ , YZ); 
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end; 



(*full spectrum*) 

X := round( 2*( I/MXDIV2 * 256 + 30)); 

Y := 300 - round(sqrt(sqr(XRA[ I] ) + sqr(XI a[ I] ) ) ) ; 
line(X,300,X,Y); 

end; 

(*graph labels*) 

setcolor( lightgray) ; 
outtextXY( 60 , 320 , ' 0 ' ) ; 
outtextXY( 331 , 320 , ' kl ' ); 
outtextXY( 391 ,320 , ' k2 1 ); 
outtextXY(565 , 320 , ' kx/2 ' ); 

outtextXYC 150,330 . ' frequency >' ); 

outtextXYC 20, 250, amplitude' ); 
readln; 

outtextXY( 300 , 330 Press return to continue. '); 

readln; 

CloseGraph; 
end; (*GraphData*) 

(* 

(* 

begin (*main body*) 

INVERSE := false; (*sets forward FFT* ) 

ERROR := 0; 

NUMPTS := kx; 

kl := round(kx * 67. 0 / 256.0 + 1); 
k2 := round(kx * 83.0 / 256.0); 

repeat 

writeln( ' Enter 1 to sample data'); 
writelnC' 2 to read data from an ASCII file'); 

readln( ANSWER); 
until ANSWER in [ ' 1 ' , ' 2 ' ] ; 

if ANSWER = V then 

(*AcquireData samples input analog signal*) 

AcquireData 
else 

GetData; 

(*GetData reads a file of discrete sample produced in the transmitter. 
GetData is used for testing Graphics and Decoding procedures without 
the need for external hardware and Data aquisition board 
initialization*) 

writeln( ' Computing FFT'); 

ComplexFFT ( NUMPTS , INVERSE , XREAL , XI MAG , ERROR ) ; 

Dif fDecode; 

(* writelnC 'Error = ', ERROR,' hit the enter key '); readln; *) 
GraphDataC XREAL, XIMAG); 
end. 



■*) 

■*) 
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APPENDIX J. RECMES 



program RECEIVER; 

(-'Acquires the signal. Stores it in a memory buffer. Differential decodes 
between tones. Maximum number of bauds are received. The number of 
bauds processed is a user input*) 

uses Graph, Crt, tp4dl6; 



(*$I-*) 

(*$R-*) 

const MAX.BUFFER = 65500; 
type 

TNvector = array-0. . 4096- of real; (-'TYPE for real and imaginary data 

for FFT routing*) 

TNvectorPtr = ATNvector; (*Pointer for FFT data array which 

allows dynamic allocation of memory*) 



var 

INVERSE .-boolean; 

XREAL, XI MAG : TNvectorPtr; 
ERROR, TEMPBYTE: byte; 

J, 

kl,k2,kx, ANSWER, 

ERR.CODE, BAUDCOUNT, 



SYMBOLCOUNT, 

NUMBAUDS , 

MAXNUMBAUDS : integer; 

MAGNITUDE, 

PHASE : real; 

DATAVECTOR : a integer; 

DMAPOINTER -.pointer; 

TESTFILE, 

OUTFILE : TEXT; 

(*$I FFT87B2. INC*) 

(*$I COMPFFT. INC*) 

(* *) 

(* *) 



procedure PacketSetUp; 



begin 

repeat 

clrscr; 

if kx < 0 then writeln('TRY AGAIN 1 ); 
writeln( ' Enter baud size '); 
readln( ANSWER) ; 
case ANSWER of 
256: kx: = 256; 

512: kx: = 512; 

1024: kx: =1024; 
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2048: kx: =2048; 

4096: kx: =4096; 
end; (*case*) 
if kx = 0 then kx := -1; 
until kx > 0; 

MAXNUMBAUDS := trunc( (MAX_BUFFER/2)/kx); 



repeat 

writeln; 

writelnC ' Enter number of ',kx,' bauds to process. ', 
MAXNUMBAUDS,' is the maximum.'); 
r ead ln( NUMB AUDS ) ; 

until NUMB AUDS in [ 1. . MAXNUMBAUDS] ; 



kl := round(kx * 67.0 / 256.0 + 1); 
k2 := round(kx * 83.0 / 256.0); 
end; (*Packet SetUp*) 

(* 

(* 

procedure AcquireData; 

(*AcquireData initializes Metrobyte DASH-16F data acquisition board, 
using TT00LS procedure D16_int and D16_ainm. Data transfer is 
controlled by the DMA controller and initialized by D16_ainm and 
disabled by D16_dma_int_disable. TT00LS procedures are external 
procedures included by 'uses' tp4dl6. *) 



-*) 

-*) 



RATE: real; 

I , CNT_NUM , MODE, CYCLE, TRIGGER, 

BASE_ADR , INT_LEVEL, DMA_LEVEL, 

BOARD NUM, CHANL0, 

0P_TYPE , STATUS, NEXT_CNT, ERR_C0DE_S : integer; 

begin 

BOARD NUM := 0; INT LEVEL := 7; DMA_ LEVEL := 1; 

BASE_ADR := $300; 

D16_init( B0ARD_NUM , BASE_ADR , INT_LE VEL , DMA_LEVEL ,ERR_C0DE ) ; 

CHANL0 := 0; 

CYCLE: =0; (*0-one sweep of the DMA 1-autoinitialize*) 

TRIGGER: =0; (*0 - external 1 - internal*) 

CNT_NUM: =trunc(MAX_BUFFER / 2); (*# of samples*) 

RATE := 10000. 0; (*used for internal trigger*) 

MODE := 2; (*DMA mode*) 
writelnC 'Ready to acquire'); 

D16_ainm( B0ARD_NUM , CHANL0 , MODE , CYCLE , TRIGGER , CNT_NUM , 

RATE , DATAVECTORa , ERR_C0DE ) ; 

STATUS := 11; 

(*status indicates the progress of acquisition. When all 
samples have been acquired status=0*) 
repeat 

D 16_dma_int_s t atus ( B0ARD_NUM , 0P_TYPE , STATUS , NEXT_CNT , ERR_C0DE_S ) ; 
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until STATUS = 0; 
writeln( ' Data received'); 

if ERR_CODE <> 0 then 

D 1 6_p r int_er r or ( ERR_C0DE ) ; 

D 1 6_dma_ int_d is ab 1 e ( B0ARD_NUM , ERR_C0DE ) ; 



end; (*Acquire*) 

(* - *) 

(* *) 

procedure ConvertData; 

(*ConvertData seperates channel and acquired data. CHAN_DATA is not used. 
Acquired data is stored in XREAL.*) 

var 

AD_DATA: arrayf 0. . 4096] of integer; 

I , CHAN_DATA , ERR_C0DE , 

SEGMENTPART , OFFSETPART : integer; 

NEWDATA VECTOR : a integer; 

TEMPPOINTER : pointer; 

begin 

writeln( ' Processing baud ' ,BAUDC0UNT); 
fillchar(XREALA ,sizeof (XREALa) , 0); 
fillchar(XIMAGA , sizeof (XIMAGa) ,0); 

SEGMENTPART: =seg(DATAVECTORA); 

OFFSETPART: =of s(DATAVECTORa ) + 2 * kx * (BAUDCOUNT - 1); 

TEMPPOINTER: =pt r ( SEGMENTPART , OFFSETPART) ; 

NEWDATA VECTOR : = TEMPPOINTER; 

dl6_convert_data( 2047 ,kx , NEWDATA VECTOR a ,AD_DATA[ 0] , 

CHAN_DATA , 0 , ERR_C0DE ) ; 



for I: = 0 to (kx - 1) do 
begin 

XREALa[I] : = AD_DATA[ I] /5; 
end; 

end; (*ConvertData*) 

(* - -*) 

procedure DiffDecode; 



(*Dif fDecode differentially decodes complex frequency domain arrays XREAL 
and XIMAG. Four decoded symbols are recombined into a byte and 
transferred to file BYTESOUT.DAT.*) 



var 

I : integer; 

TEMPRE AL , TEMP IMAG .-real; 
BITS : byte; 

TEMPCHAR t char; 



begin 

for I : = kl to (k2 - 1) do 
begin 

(*Complex multiply two adjacent tones, I and the complex conjugate of 1+1. 
This will give the phase difference between the two tones. The answer 
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is in rectangular notation*) 

TEMPREAL: =XREALa[ I] * XREALa[ 1+1] + XIMAGa[ I] * XIMAGa[ 1+1] ; 
TEMPIMAG: =XREALa[ I] * XIMAGa[ 1+1] - XREALa[ 1+1] * XIMAGa[ I] ; 

(*Complex multiply (TEMPREAL + j TEMPIMAGE) and (1+j). This rotates the 
differential vector pi/4 radians. XREAL and XIMAG are used to store 
the results. This eliminate the original data*) 

XREALa[ I] :=( TEMPREAL - TEMPIMAG)/80; 

XIMAGa[ I] :=( TEMPREAL + TEMPIMAG)/80; 



(*decode*) 

if (XREALa[I] >= 0) and (XIMAGa[ I] > 0) then 
BITS: =$00; 

if (XREALa[I] < 0) and (XIMAGa[ I] > 0) then 
BITS: =$01; 

if (XREALa[I] < 0) and (XIMAGa[ I] <= 0) then 
BITS: =$03; 

if (XREALa[I] >= 0) and (XIMAGa[ I] <= 0) then 
BITS: =$02; 

SYMB0LC0UNT := SYMB0LC0UNT + 1; 



(*fill TEMPBYTE with four symbols*) 
if f rac( SYMB0LC0UNT / 4) = 0.25 then 
TEMPBYTE := (BITS shl 6); 
if f rac( SYMB0LC0UNT / 4) = 0. 5 then 
TEMPBYTE := (BITS shl 4) or TEMPBYTE; 
if f rac( SYMBOLCOUNT / 4) = 0.75 then 
TEMPBYTE := (BITS shl 2) or TEMPBYTE; 
if (frac( SYMBOLCOUNT / 4) = 0.0) then 
begin 

TEMPBYTE := BITS or TEMPBYTE; 

TEMP CHAR := chr( TEMPBYTE); 

write( 0UTFILE ,TEMPCHAR) ; 

TEMPBYTE: =0; 
end; (*if frac*) 
end; (*for I*) 



XREALa [ k2] : =1. 0; 
XIMAGa[ k2] : =1. 0; 



end; (*Dif fDecode*) 

(* *) 

procedure Showmessage; 

(*Showmessage reads in decoded message*) 



NEXTCHAR: char; 

begin 

write In; 

writeln('The message transmitted is..'); 
assign(0UTFILE, ’MESSAGE. DAT’ ); 
reset(OUTFILE); 
while not EOF(OUTFILE) do 
begin 

while not EOLN(OUTFILE) do 
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begin 

read(OUTFILE,NEXTCHAR); 
write(NEXTCHAR); 
end; ( *wh i 1 e not EOLN*) 
readln( OUTFILE); 
writeln; 

end; (*while not EOF*) 
close( OUTFILE); 
end; (*Showmessage*) 

(* *) 

begin (*raain body*) 

GetDMABuf fer( MAX_BUFFER , DMAPOINTER , ERR_CODE ) ; 

DATAVECTOR := DMAPOINTER; (*This statement assigns a generic pointer 
to a variable of a specific pointer type, 
i. e. Ainteger, so that the pointer can be 
passed to the d!6_ainm routine.*) 



new(XREAL); 

new(XIMAG); 



INVERSE := false; (*sets forward FFT*) 

ERROR := 0; 
kx: =0; 



PacketSetUp; 

(*received message output file*) 
ass ign( OUTFILE, 'MESSAGE. DAT' ); 
rewrite( OUTFILE); 

(*RECSTAT file shows the computed real and imaginary values of the 
received signal and its decoded representation*) 
ass ign( TESTFILE , ' RECSTAT. DAT' ) ; 
rewrite (TESTFILE); 

(* writeln(TESTFILE, ' Baudlength is',kx:5);*) 

SYMB0LC0UNT: =0; 

TEMP BYTE: =0; 

AcquireData; (*AcquireData samples input analog signal*) 

for BAUDCOUNT : = 1 to NUMBAUDS do 
begin 

ConvertData; 

ComplexFFT (kx, INVERSE, XREAL,XIMAG. ERROR); 

(* writeln(TESTFILE, 'Baud': 11 'Tone : 5 'Real' : 10, 

' Imaginary' : 10, 'Mag' : 10, 'Phase' : 8); *) 
for J:=kl to k2 do 
begin 

MAGNITUDE: =20*ln( sqrt( sqr(XREALA[ J] )+ 
sqr(XIMAGA[ J] ))) / ln( 10); 

if (XREALa[J] >= 0) then 
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PHASE: =arctan(XIMAGA[J] / XREALa[J]); 
if (XREALa[J] < 0) and (XIMAGa[J] > 0) then 
PHASE: =Pi + arctan(XIMAGA[ J] / XREALa[J]); 
if (XREALa[J] < 0) and (XIMAGa{J] <= 0) then 
PHASE: =arctan(XIMAGA[ J] / XREALa[J]) - Pi; 

writeln(TESTFILE , ' R 1 ,BAUDCOUNT: 4, J: 5, 

XREALa[ J] : 10: 4,XIMAGa[ J] : 10: 4, 
MAGNITUDE: 10: 4, PHASE: 8: 4); 

end; 

Dif fDecode; 

for J:=kl to k2 do 
begin 

MAGNITUDE: =sqrt(sqr(XREALA[J] )+sqr(XIMAGA[ J] )); 
WritelnCTESTFILE/D' , BAUDCOUNT: 4 , J: 5 , 

XREALa[ J] /MAGNITUDE: 10: 4, 
XIMAGa[ J] /MAGNITUDE: 10: 4); 

end; 

end; 

writeln(OUTFILE , chr(TEMPBYTE) ) ; 

close( OUTFILE) ; 
close(TESTFILE); 
dispose(XREAL); 
dispose(XIMAG); 

F r eeDMABuf f e r ( MAX_BUFFER , DMAPOINTER , ERR_C0DE ) ; 

Showmessage; 

(* writeln( 'Error = 1 , ERROR, ' hit the enter key ' ); readln; *) 
end. 
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APPENDIX K. RESPONSE 



program RESPONSE; 

(*RESPONSE outputs the channels phase response, by subtracting 
the transmitted phase from the received phase for each tone. 
Prior to running, set NUMTONES and K1MINUS1 for the proper 
baud parameters then compile and run the program*) 



MAG, PHASE, 
XREAL,XIMAG 
K1MINUS1 , 
NUMTONES, K, 

I ,BAUD,Kx,Kr 
XMITPHASE, 
XR,XI 
CODE 

XMITDATA, 
RECDATA , 
CHANNELDATA 



begin 

ass ign( XMITDATA, 'xmitdatl. 024' ); 
reset (XMITDATA); 

as sign( RECDATA, ' recstat. dat' ); 
reset( RECDATA); 

as sign( CHANNELDATA, ' resp 1024. dat' ); 
rewrite( CHANNELDATA); 

(*NUMT0NES and K1MINUS1 should be set for the same baud parameters 
as xmitdat*. *** and recstat. dat*) 

NUMTONES: =64; 

K1MINUS1: =268; 

(*get transmitted phase and received phase*) 
while not EOF (XMITDATA) do 
begin 

for I: =1 to NUMTONES do 
begin 

read ln( XMITDATA, BAUD, Kx,XREAL[ I] ,XIMAG[ i] ); 
readln( RECDATA, CODE, BAUD, Kr,XR, XI, MAG[ I] ,PHASE[ I] ); 
end; 

(*ensure tones match*) 
if Kx <> Kr then 
begin 

writeln('Kx <> Kr'); 
readln; 
end; 

(*convert transmitted symbol into phase*) 



: array[ 1. . 256] of real; 

: integer; 

: real; 

: char; 

: text; 
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for I: =1 to NUMTONES do 
begin 

if (XREAL[ I] >= 0.0) and (XIMAG[ I] >= 0.0) then 
XMITPHASE: = pi/4; 

if (XREAL[ I] < 0.0) and (XIMAG[i]>= 0.0) then 
XMITPHASE: = 3*pi/4; 

if (XREAL[ I] >= 0.0) and (XIMAG[ I] < 0.0) then 
XMITPHASE: = -pi/4; 

if (XREAL[ I] < 0.0) and (XIMAG[I]< 0.0) then 
XMITPHASE: = -3*pi/4; 

(^subtract transmitted phase from received phase*) 

PHASE [ I] : =PHASE[ I] - XMITPHASE; 

(*phase differences may jump be 2*pi radian. This helps 
but does not work in every case. Output file may have 
to be edited to 2*pi jumps*) 

0 if PHASE[ I] > 2*pi then 

PHASE[ i] : =PHASE[ I] - 2*pi; 
if PHASE[ i] < -2 then 

PHASE[ i] : =PHASE[ i] + 2*pi; 

(*actual tone numbers*) 

K:=I + K1MINUS1; 

writeln(CHANNELDATA,BAUD: 4,K: 4,MAG[ I] : 10: 4,PHASE[ I] : 10: 4); 
readln(RECDATA); 
end; 

end; 

close(XMITDATA); 

close(RECDATA); 

close(CHANNELDATA); 

end. 
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APPENDIX L. STATISTICS 



program Statistics; 

(*to run this program the transmitter program Xmitmes must be run using 
MESSAGE.DAT. The file produced, XMITDAT.DAT, must be transferred to the 
receiver computer. Then transfer the data between computers using 
Xmitmes and Receive. Receive produces the file RECSTAT.DAT. This 
program uses XMITDAT. DAT and RECSTAT. DAT to compute the mean and 
variance for each quadrant, and the SNR for each quadrant, baud and 
the overall SNR for the packet.*) 

uses CRT; 

type 

STATARRAY = array[ 1. . 125] of real; 



Kx,K,NUMBAUDS , 

I, J, BAUD, 

TONEX, TONER, 

Q1,Q2,Q3,Q4, 

BITERR0R1 ,BITERR0R2 , BITE RR0R3, BITE RR0R4 : integer; 

TEMPR, TEMPI, 

SNRINDB , SNRIN , 

SNROUTDB , SNROUT, 

MEANR , MEANI , VARR , VARI : r ea 1 ; 

CODE : char; 

XD,RD,SD : text; 

XR,XI ,RR,RI : array[ 1. . 16,1. . 128] of real; 

BTM , BTV , BTQ : array[ 1. . 16] of real; 

SRI , SR2 ,SR3 , SR4 , 

SI1,SI2,SI3,SI4 : STATARRAY; 

STATMATRIC : array[ 1. . 16 , 1. . 4, 1. . 5] of real; 



(* *) 

procedure Stat(NPTS: integer; 

X: STATARRAY; 

var XMEAN , XVAR :real); 

(*input NPTS: number of points to be used to compute the mean 
X: the array of data 

output XMEAN: the mean of the data in array X 

XVAR: the variance of the data in array X*) 



SUM: real; 

N : integer; 
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begin 

(*Compute mean*) 

SUM: = 0. 0; 

for N:= 1 to NPTS do 
SUM: = SUM + X[ N] ; 

if NPTS > 0 then 

XMEAN: = SUM / NPTS 
else 

XMEAN: =1. 0e-10; 

(*Compute variance*) 

SUM: = 0. 0; 

for N: =1 to NPTS do 

SUM: = SUM + sqr(X[N] - XMEAN); 

if NPTS > 1 then 

XVAR: = SUM / (NPTS - 1) 
else 

XVAR: =1. Oe- 10; 
end; (*Stat*) 

(*- --*) 

begin (*main body*) 
repeat 

writeC'Enter baud size. '); 
readln(Kx); 

until (Kx mod 256 = 0) and (Kx >= 256) and (Kx <= 4096); 

writeC'Enter input SNR in DB ' ); 
readlnC SNRINDB) ; 

SNRIN := exp( SNRINDB / 10 * ln(10)); 

K := trunc(Kx / 16); 

as s ign( XD , ' XMITDAT. DAT ' ) ; 
reset(XD); 

as s ign( RD , ' RECSTAT. DAT' ); 

reset(RD); 

case Kx of 

256: begin assign(SD, ' STAT256. DAT' ); NUMBAUDS 
512: begin assign( SD , ' STAT512. DAT' ) ; NUMBAUDS 
1024: begin assign( SD, ' STAT1024. DAT' ); NUMBAUDS 
2048: begin assign( SD , ' STAT2048. DAT' ); NUMBAUDS 
4096: begin ass ign( SD , ' STAT4096. DAT ' ) ; NUMBAUDS 
end; (*case Kx*) 
rewrite(SD); 

(*read in transmitted and received data*) 
for I: = 1 to NUMBAUDS do 
begin 

for J: = 1 to K do 
readln(RD); 



— ±o; ena; 
= 8; end; 
= 4; end; 
= 2; end; 
= 1; end; 
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for J: = 1 to K do 
begin 

readln(XD, BAUD, TONEX, XR[I,J], XI[I,J]); 
readln(RD, CODE, BAUD, TONER, RR[I,J], RI[I,J]); 
if TONEX <> TONER then 
begin 
clrscr; 

writelnCXMITDAT.DAT and RECSTAT.DAT are in error. '); 
Halt; 
end; 

end; (*for J*) 
end; (*for I*) 
close(XD); 
close(RD); 



(*decode transmitted data*) 
for I: = 1 to NUMBAUDS do 
for J;= 1 to (K - 1) do 
begin 

TEMPR: = XR[ I , J] 
TEMPI: = XR[ I , J] 
XR[ I , J] : = TEMPR 
XI[I,J] := TEMPR 
end; (*for J*) 



* XR[I,J+1] +XI[I,J] * XI [ I , J+l] ; 

* XI[I,J+1] - XR[ I , J+l] * XI [ I , J] ; 
- TEMPI; 

+ TEMPI; 



(*store statistics for each quadrant and baud compute statistics*) 
BITERR0R1: =0; BITERR0R2: =0; BITERR0R3: =0; BITERR0R4: =0; 
for I: = 1 to NUMBAUDS do 
begin 

Ql: =0; Q2: =0; Q3: =0; Q4: =0; 
for J: = 1 to (K - 1) do 
begin 

if (XR[ I , J] >= 0.0) and (XI[I,J] > 0.0) then 
begin 

Ql: =Q1 + 1; 

SR1[ Ql] := RR[ I , J] ; 

SI 1[ Ql] := RI[ I , J] ; 
if RR[ I , J] < 0. 0 then 

BITERR0R1: = B I TERROR 1 + 1; 
if RI[ I , J] <= 0. 0 then 

BITERR0R1: = BITERROR 1 + 1; 

end; 

if ( XR[ I , J] < 0.0) and (XI[ I , J] > 0.0) then 
begin 

Q2:=Q2 + 1; 

SR2[Q2] := RR[ I , J] ; 

SI2[Q2] := RI[ I , J] ; 
if RR[ I , J] >= 0. 0 then 

BITERR0R2: = BITERR0R2 + 1; 
if RI[ I , J] <= 0. 0 then 

BITERR0R2: = BITERR0R2 + 1; 

end; 

if (XR[ I , J] < 0.0) and (XI[I,J] <= 0.0) then 
begin 

Q3:=Q3 + 1; 

SR3[Q3] := RR[ I , J] ; 
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SI3[Q3] := RI[ I , J] ; 
if RR[ I , J] >= 0. 0 then 

BITERR0R3: = B I TERROR 3 + 
if RI[ I , J] > 0. 0 then 

BITERR0R3: = BITERR0R3 + 

end; 

if (XR[ I , J] >= 0.0) and (XI[ I , J] 
begin 

Q4:=Q4 + 1; 

SR4[ Q4] : = RR[ I , J] ; 

SI4[ Q4] : = RI[ I , J] ; 
if RR[ I , J] < 0. 0 then 

BITERR0R4: = BITERR0R4 + 
if RI[ I , J] >0.0 then 

BITERR0R4: = B I TERROR 4 + 

end; 

end;(*for J*) 



1; 

1 ; 

<= 0. 0) then 



1 ; 

1 ; 



(*STATMATRIC is a 3 
number. The second 
for each quadrant. 
STATMATRIC[ I , J , 1] = 
STATMATRIC[ I , J,2] = 
STATMATRIC[ I , J, 3] = 
STATMATRIC[ I , J , 4] = 
STATMATRIC[ I , J ,5] = 



dimensional array. The first index is the baud 
is the quadrant. The third is the statistics 

; the mean of the real parts 
: the variance of the real parts 
: the mean of the imaginary parts 
: the variance of the imaginary parts 
: the number of points transmitted in a given 
quadrant*) 



BTM[ I] : =0. 0; 

BTV[ I] : =0. 0; 

BTQ[ I] : =0. 0; 
if Q1 > 1 then 
begin 

Stat( Ql, SRI, MEANR,VARR); 

Stat( Q1 , SI1 ,MEANI , VARI); 

STATMATRIC[ 1,1,1]: =MEANR; 

STATMATRIC[ 1,1,2]: =VARR; 

STATMATRIC[ 1,1,3]: =MEANI; 

STATMATRIC[ 1,1,4]: =VARI; 

BTM[ I] : =BTM[ I] + Q 1 * ( ab s e ME ANR ) +ab s e ME AN I ) ) ; 
BTV[ I] : =BTV[ I] + eQl"l)*(VARR+VARI); 

BTQ[ I] : =BTQ[ I] + 2*Q1; 
end; 

STATMATRIC[ 1,1,5] : =Q1; 



if Q2 > 1 then 
begin 

StateQ2,SR2,MEANR,VARR); 

Stat( Q2,SI2,MEANI, VARI); 

STATMATRIC[ 1,2,1]: =MEANR; 

STATMATRIC[ 1,2,2] : =VARR; 

STATMATRIC[ 1,2,3]: =MEANI; 

STATMATRIC[ 1,2,4] :=VARI; 

BTM[ I] : =BTM[ I] + Q2*eabseMEANR)+abseMEANI) ); 
BTV[ I] : =BTV[ I] + eQ2-l)*(VARR+VARI); 

BTQ[ I] : =BTQ[ I] + 2*Q2; 
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end; 

STATMATRIC[ 1,2,5] :=Q2; 

if Q3 > 1 then 
begin 

Stat ( Q3 , SR3 , MEANR , VARR) ; 

Stat(Q3 , SI3 ,MEANI ,VARI); 

STATMATRIC] 1,3,1]: =MEANR; 

STATMATRIC] 1,3,2]: =VARR; 

STATMATRIC[ 1,3,3] : =MEANI; 

STATMATRIC] 1,3,4]: =VARI; 

BTM[ I] : =BTM[ I] + Q3*( abs ( MEANR )+abs( MEAN I ) ) ; 
BTV[ I] : =BTV[ I] + (Q3-l)*( VARR+VARI); 

BTQ[ I] : =BTQ[ I] + 2*Q3; 
end; 

STATMATRIC[ 1,3,5] : =Q3; 

if Q4 > 1 then 
begin 

Stat( Q4 , SR4 , MEANR , VARR) ; 

Stat(Q4,SI4,MEANI ,VARI); 

STATMATRIC] 1,4,1]: =MEANR; 

STATMATRIC] 1,4,2] : =VARR; 

STATMATRIC[ 1,4,3] : =MEANI; 

STATMATRIC] 1,4,4] :=VARI; 

BTM[ I] : =BTM[ I] + Q4*( abs ( MEANR ) +abs ( ME AN I ) ) ; 
BTV[ I] : =BTV[ I] + (Q4-l)*( VARR+VARI); 

BTQ[ I] : =BTQ[ I] + 2*Q4; 
end; 

STATMATRIC] 1,4,5] : =Q4; 
end; (*for I*) 



(*output statistics to a file*) 

write ln(SD, '*** SNRIN= ' ,SNRIN: 10: 7, ' or ' , SNRINDB: 10: 7, ' DB ***'); 
writeln(SD); 

writeln(SD,' Baud length is ' ,Kx); 



SR2[ 1] : =0. 0; 

SR2[ 2] : =0. 0; 

SR2[ 3] : =0. 0; 

for I: = 1 to NUMBAUDS do 
begin 

writeln( SD); write ln( SD); write ln( SD) ; write ln( SD) ; 
writeln( SD, 1 DFT statistics for baud #; ' , I: 4); writeln( SD); 

for J: = 1 to 4 do 
begin 

(*real output data*) 

writeln(SD); write ln(SD); 

writeln( SD, ' Given the phase is in quadrant ' , J); 
writeln( SD); 

if STATMATRIC] I , J,5] > 1 then 
begin 

writeln(SD, 'Number of points = T : 35 , 
STATMATRIC[ I , J,5] : 10: 5); 
writeln(SD, 'Mean of real parts = ' : 35 , 
STATMATRIC[ I , J, 1] : 10: 5); 
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writeln(SD, 'Variance of real parts = ’ : 35 , 
STATMATRIC[ I , J, 2] : 10: 5); 

SNROUT: =sqr( STATMATRIC[ I , J, 1] ) / STATMATRIC[ I , J,2] ; 
SNROUTDB: = 10.0 * ln(SNROUT) / ln(10.0); 
writeln(SD, ' SNROUT of real part = ' : 35 , 

SNROUT: 10: 5,' = ' , SNROUTDB: 10: 5 , ' DB'); 
writeln(SD, ' SNROUT/SNRIN of real part = ' : 35 , 
(SNROUT / SNRIN) : 10: 5); 

(^imaginary output data*) 
writeln(SD); 

writeln(SD, 'Mean of imag parts = ' : 35 , 

STATMATRIC[ I , J, 3] : 10:5); 
writeln( SD, ' Variance of imag parts = ' : 35 , 
STATMATRIC[ I , J,4] : 10: 5); 

SNROUT: = sqr(STATMATRIC[ I , J, 3] ) / STATMATRIC[ I , J,4] ; 
SNROUTDB: = 10 * ln(SNROUT) / ln(10.0); 
writeln( SD, ' SNROUT of imag part = ' : 35 , 

SNROUT: 10: 5, ' = ', SNROUTDB: 10: 5 , ' DB’); 
writeln(SD,' SNROUT/SNRIN of imag part = ' : 35 , 
(SNROUT / SNRIN): 10: 5); 

end (*if Statmatric*) 

else 

begin 

writeln(SD, ' Quadrant has less then 2 statistic 
points '); 

writeln(SD, ' and is not used in the overall statics.'); 
end; (*else*) 
end; (*for J*) 

write ln( SD); write ln( SD); write ln( SD); 

writeln( SD, ' Overall (Real + Imag) statistics for Baud #: ',1:4); 
SNROUT: = sqr(BTM[ I] /BTQ[ I] ) / (BTV[ I] /(BTQ[ I] -1)); 

SNROUTDB: = 10 * ln(SNROUT) / ln(10.0); 

writeln(SD, 'Overall Baud SNROUT = ' : 35 , SNROUT: 10: 5 , ' = ', 
SNROUTDB: 10: 5, ’ DB'); 

writeln(SD, 'Overall baud gain (SNROUT / SNRIN) = ' : 35 , 
(SNROUT/SNRIN): 10: 5); 

(*store baud means and variances*) 

SR2[ 1] := SR2[ 1] + BTM[ I] ; 

SR2[ 2] : = SR2[ 2] + BTV[ I] ; 

SR2[ 3] : = SR2[ 3] + BTQ[ I] ; 

end; (*for I*) 

writeln( SD); write ln( SD); writ eln(SD); write ln( SD) ; 

writeln( SD , ' *** TOTAL OVER ALL ' , NUMBAUDS , ' BAUDS ABOVE BAUD ', 

Kx , ' *** ' ) ; 

writeln( SD) ; write ln( SD) ; 

SNROUT: =sqr(SR2[l]/SR2[ 3] ) / (SR2[2] / ( SR2[ 3] - 1)); 

SNROUTDB: =10 * ln(SNROUT) / ln(10.0); 

writeln(SD, 'Overall SNROUT = ' : 35 , SNROUT: 10: 5 , ' = ', 

SNROUTDB: 10: 5, ’ DB’ ); 
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writeln( 'Overall SNROUT = ' : 35 , SNROUT: 10: 5 , ' = ', 

SNROUTDB: 10: 5, ' DB'); 

writeln(SD, 'Overall gain (SNROUT / SNRIN) = ' : 35 , 

(SNROUT/SNRIN): 10: 5); 

writeln( SD) ; writeln(SD); 

writeln( SD, 'Total possible statistical points ' , (Kx/8-2)*NUMBAUDS: 3: 0, 
'. Actual number of statistical points ' , SR2[ 3] : 3: 0) ; 



writeln( SD); 

writeln(SD, 'BIT ERRORS': 35); 
write ln( SD, ' Quadrant 1 = ' : 30 , BITERR0R1); 
writeln( SD, ' 2 = ' : 30,BITERROR2); 
writeln( SD, ' 3 = ' : 30,BITERROR3); 
writeln( SD, ' 4 = ' : 30 ,BITERR0R4); 

B I TERROR 1: = B I TERROR 1 + BITERR0R2 + BITERROR3 + 
writeln(SD, 'Total Bit errors = ' : 35 , BITERR0R1); 
writeln( 'Total Bit errors = ':30, BITERR0R1); 



BITERR0R4; 



close(SD); 
end. (*main body*) 
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